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ABSTRACT 

We use recent progress in simulating the production of magnetohydrodynamic jets 
around black holes to derive the cosmic spin history of the most massive black holes. 
I Our work focusses on black holes with masses ^ 10® M©. Under the assumption that 

HH , the efficiency of jet production is a function of spin a, as given by the simulations, we 

. • can approximately reproduce the observed 'radio loudness' of quasars and the local 

' radio luminosity function. 

Using the X-ray luminosity function and the local mass function of supermassive 
Q ' black holes, SMBHs, we can reproduce the individual radio luminosity functions of 

i—{ , radio sources showing high- and low-excitation narrow emission lines. We find that 

■ the data favour spin distributions that are bimodal, with one component around spin 
' zero and the other close to maximal spin. The 'typical' spin is therefore really the 

expectation value, lying between the two peaks. In the low-excitation galaxies, the two 
components have similar amplitudes, meaning approximately half of the sources have 
^ I very high spins, and the other have very low spins. For the high-excitation galaxies, 

. the amplitude of the high-spin peak is typically much smaller than that of the low- 

PsJ ' spin peak, so that most of the sources have low spins. However, a small population of 

CS| ! near maximally spinning high-accretion rate objects is inferred. A bimodality should 

O^l ' be seen in the radio loudness of quasars, although there are a variety of physical and 

' selection effects that may obfuscate this feature. 

, We predict that the low-excitation galaxies are dominated by SMBHs with masses 

^10® Mq, down to radio luminosity densities '^10^^ W Hz~^ sr~^ at 1.4 GHz. Under 
reasonable assumptions, our model is also able to predict the radio luminosity function 
at z=l, and predicts it to be dominated by radio sources with high-excitation narrow 
^ emission lines above luminosity densities J^IO^^ W Hz^^ sr^^ at 1.4 GHz, and this is 

^ I in full agreement with the observations. 

■ From our parametrisation of the spin distributions of the high- and low-accretion 

rate SMBHs, we derive an estimate of the spin history of SMBHs, which shows a 
weak evolution between z=l and 0. A larger fraction of low-redshift SMBHs have 
high spins compared to high-redshift SMBHs. Using the best fitting jet efficiencies 
there is marginal evidence for evolution in spin: the mean spin increases slightly from 
(0)^^0.25 at z=l to (a)~0.35 at z—0, and the fraction of SMBHs with increases 
from 0.16±0.03 at z=l to 0.24±0.09 at z^O. Our inferred spin history of SMBHs is 
in excellent agreement with constraints from the mean radiative efficiency of quasars, 
as well as the results from recent simulations of growing SMBHs. We discuss the 
implications in terms of accretion and SMBH mergers. We also discuss other work 
related to the spin of SMBHs as well as work discussing the spin of galactic black 
holes and their jet powers. 

Key words: galaxies : active galaxies : jets-galaxies: nuclei - quasars: general -black 
hole physics - cosmology: miscellaneous 
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1 INTRODUCTION 

Astrophysical black holes are described by two parameters, 
mass and spin, which define the structure of space-time 
around the black hole. Distant observers can infer black hole 
properties from the dynamics of material orbiting the black 
hole or more indirectly when matter is being accreted, from 
the radiation and jets produced. 

Active galactic nuclei (AGN) are believed to be pow- 
ered bi;_accretion onto supermassi ve black holes (SMBHs, 
e.g. lLvnden-Belllll969l : iReed Il984l ). Matter accreted onto 
the SMBH radiates via thermal and non-thermal processes, 
dominating the observed spectral energy distribution (SED) 
bet ween X-ray e nergies and far-infrared wavelengths (see 
e.g. lKrolik|[l999l ). Multiwavelength observations can there- 
fore be used to measure the bolometric luminosity of the 
AGN, although in practice reasonable estimates of this 
qua ntity can be obtained from observation s in a single band 
(e.g. lElvis et al.l[l993 : iHopkins et al]l2007l l . 

Accretion onto SMBHs can also produce jets. The bulk 
of the jet power is used to fill and expand a cocoon of 
low-density relativistic magnetised plas ma that pushes back 
the intergalactic medium (IGM, see e.g. iLongair et aPl 19731 : 
IScheuei<1973 : lBegel man et al]|l984l ). These lobes are sources 
of synchrotron radiation emitted at radio wavelengths. How- 
ever, the power is radiated at a much lower rate than is 
supplied to the lobes via the jets, at least while the jets are 
active, so that energy is stored in the lobes. Making assump- 
tions about the contents of the lobes (charged particles and 
magnetic fields), and how the energy density is distributed 
between them, as well as assumptions on the typical jet ad- 
vance speeds and IGM density profiles, the jet power can 
be estimated fr om the low-frequency monochromatic radio 
luminosity (e.g. iMiller et al ] |l993l : IWillott et al.l[l999l ). 

Observations of AGN require an extra 'hidden variable' 
to explain the variations in radio luminosity, reflecting vari- 
ations in jet power. For a given radiative power, AGN can 
show radio lumi nosities that vary by mor e than two orders of 
magnitude (e.g. iKellermann et aDll989l: filler et al.lll993l : 



ISikora et alll2007l : iMcNamara et al.ll201]} ). This huge scat- 
ter in 'radio-loudness' is not dependent on black hole mass 
or accretion rate (e.g. ISikora et al]l2007l ). and the X-ray to 
far-infrared SEDs of radio-loud and radio-quiet AG N do not 
show significant differences (e.g. lElvis et al.l|l994 ). 

Gonversely, observational studies of AGN selected at 
radio frequencies have, for a long time, realised that the 
AGN showed two types of optical spectra. Some show high- 
excitation narrow emission lines (known as 'high-excitation 
galaxies' or HEGs) indicative of high accretion rates, while 
the others show only low-excitation narrow lines (LEGs) 
suggesting low ac c retion rates (e.g. Hine fc Longaiij 1 19791 : 
iLaing et al.l Il983l . Il994l : I Jackson fc Rawlingd 1 19m This 
confirms that SMBHs producing jets are observed to accrete 
both at high and low rates. 

Gomparison of the jet powers and bolometric luminosi- 
ties of AGN suggests that the conversion of accreted energy 
into jets must occur extremely efficiently in some sources, 
since the inferred jet power is comparable or so metimes even 
greater than the inferred accretion power rhc (|Punslvll2007l . 
I2OI0I : iMcNamara et al.ll201ll : iFernandes et aLlboid l. These 
results still carry considerable uncertainty in their estimates 



of jet power and accretion rate, however they do show that 
an extremely efficient mechanism must be powering the jets. 

In theory, spinning black holes can reach nominal jet 
efficiencies greater than one, since as well as drawing en- 
ergy from the accreted mass, the jet draws energy from 
the rotation of the black hole. Hence, spin is an attrac- 
tive candidate to explain t he radio loudness of AGN (e.g 
'Blandford & Znaiekl 'l9 77l: iReesI '1984: P unslv fc Goronit: 



1990; Scheuer .1992: .Wilson fc Golbert 199^ iKoide et al 
2Q0i ISikora et al.ri2007l ;i. It is attractive because, in prin- 
ciple, it is capable of high efficiencies and in addition it pro- 
vides a physical explanation for a hidden variable. This is 
usually referred to as the 'spin paradigm'. 

The jets are most likely powered by differential rota- 
tion of poloidal magnetic fields in the accretion disc or 
near the black hole. Theoretically, the spin of the black 
hole is therefore a promising candidate since frame drag- 
ging around spinning black holes would provide particularly 
strong magnetic fields an d particularly fast rotation (e.g 



Blandford fc ZnaieS Il977l : IPunslv fc Coronitil Il99(]l : iMeieil 
200ll : lKoide et al.ll2002l ). 



However, whether the jets are really tapping the rota- 
tional energy of the black hole is still a matter of debate. 
Some authors argue that the power extracted in this way 
is at most comparable to the power generated by the inner 
regi ons of the a. ccretion disc (jGhosh fc Abramowicd Il997l : 
,Livio et al.lll999l l. Others argue that in the case of geomet- 
rically thick accretion discs , powerful je ts can be be driven 
by the black hole spin (e.g. lMeieill200ll ). 

In galactic black holes (GBHs), the presence of steady 
jets occurs typically during low accretion rates (where 
by low we mean an Eddington ratcQ A<;0.01), or during 
high accretion rat es (A<;0.3, e.g. GRS 1915+105, see e.g. 
iFender et ai]|2004l '). The jets can be steady, transient or ab- 
sent at A^O.Ol, depending on the accretion state. This is 
consistent with steady jets requiring geometrically-thick ac- 
creti on flows, either AD AF s (advection-d ominated accretion 
flows iNaravan fc Yi|[l995l ) or slim discs (lAbramowicz et al.l 
Il988l l. However, these differences in the accretion flow are 
observable at X-ray energies, while as explained in the pre- 
ceding paragraphs, the radio-loudness of AGN requires an 
extra variable parameter that is not readily observable. 

Finally, we note that comparison of individual spin mea- 
surements for GBHs flnd no correlation between radio loud- 
ness and individual spin estimates dFender et al..,201 0). We 
discuss these results in detail later in this paper, in Sec- 
tion 17.21 and argue that they are not robust enough to in- 
validate the spin paradigm. 

In this work we begin by spelling out our assumptions in 
Section (2] where we also introduce the different jet efficien- 
cies, drawn from the literature, that we will use. The most 
important assumption is that the efficiency with which jets 
are produced is an increasing function of spin. In Section [3] 
we test the spin paradigm and show that under our assump- 
tions we can explain the observed range in 'radio-loudness' 
of quasars. 

In Section [4l we model the contribution to the local 



^ The Eddington rate A is defined here as A = where rh is 

the accretion rate and rrig is the accretion rate corresponding to 
the Eddington-Hmit. 



© RAS. MNRAS OOO.[Tlf30l 



Observational constraints on the spin of SMBHs 3 



radio luminosity function (LF) of HEGs and LEGs as be- 
ing composed of high-accretion rate objects (quasars or 
QSOs) and low-accretion rate ones (ADAFs). We allow 
these populations to have several different spin distributions, 
parametrised by three different models. Fitting to the data, 
and using a bayesian model selection scheme, in Section[5]we 
derive the spin distribution of SMBHs of mass M© 
in the local Universe, and make a prediction for the z=l 
radio LF. In Section |6] we present an estimate of the cos- 
mic spin history of the most massive black holes from our 
modelling of the radio LF. Section[7]compares and discusses 
our inferrence about the spins of black holes to other re- 
sults in the literature, using different methods. It includes a 
discussion about the lack of correlation between spin mea- 
surements and radio loudness found in GBHs. Finally, in 
Section [8] we discuss the implications of our result in terms 
of SMBH growth through accretion and mergers. 

We adopt a ACDM cosmology with the following pa- 
rameters: h = -H'o/(100 km s"^ Mpc"^) = 0.7; f7m = 0.3; 
f^A = 0.7. 
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Figure 1. The radiative efficiency as a function of spin. The figure 
compares the analytical model of NT73 (solid line) to the results 
from different GRMHD simul ations. The simulati on data are from 
iMcKinnev &: G ammi^ ||2004| , emp ty triangles), iBeckwith et al.l 
ll2008bl. stars'). iNoble et al.l l|2009l . black circle at a=0.9) and 
iPenna et al. empty diamonds). 



2 ASSUMPTIONS 



2.1 Radiative efficiency 



We assume that the bolometric luminosity (due to radia- 
tion), Lboi is given by: 



(1) 



where m is the accretion rate and c the speed of light. The 
te rm e is the radiative efficie ncy. In the classical disc picture 
of iNovikov fc Thornd (|l973l . hereafter NT73), e is given by 
(one minus) the binding energy per unit mass of the ISCO, 
which is itself a function of the black hole spin. 

Some simulations have shown that the radiative effi- 
ciency estimated by NT73 to slig htly underestimate the 
radiative efficienc y for low spins (jBeckwith et al.l l2008bl : 
iNoble et al.l [2OO9I ) . On the other hand, other simulations 
suggest that the radiative efh ciency is very close to that 
predicted by the NT73 model JMcKinnev fc Gammiell2004l : 
IShafee et alll2008l : iPenna et al Il201ol v 

Figure [1] compares the radiative efficiency as a function 
of spin from the NT73 model to the results of simulations 
by different groups. In some cases, the results from the sim- 
ulations have higher efficiencies than the NT73 model, in 
others they have lower efficiencies. Overall, and given the 
range of results from the simulations, the NT73 model is 
not obviously wrong. 

For our work we require analytic expressions for e(a) 
and the jet efficiency T]{a), and we will use the NT73 approx- 
imation for e(a). Figure [T] suggests this is an acceptable ap- 
proximation. This paper is concerned with variations in jet 
powers that span orders of magnitudes, and the ~2 uncer- 
tainty in e has only a small effect on our results. Indeed the 
uncertainties in the jet efBciencies discussed in Section 12.21 
are far larger than this. 

We are also neglecting the fraction of energy ad- 
vected into the hole (in excess of the energy of the ISCO). 
This is a reasonable assumption for SMBHs with Ed- 
dington rates in the range 10~^Si A ^1. Accretion rates 
with higher or lower Eddington ratios will tend to have 



a large fraction of the energy advec ted into the hole 
(jNaravan fc Yilll995l : HeloborSSIiiii). Throughout this 
paper we will only require to use the radiative efficiency for 
quasars, not for ADAFs. Quasars have typically 0.1J5A;$1.0 
(jMcLure fc Dunlod 12004 ). and hence neglecting the ad- 
vected energy is a reasonable approximation. 

2.2 Jet efficiencies 

We assume the power available for the production of jets 
can be approximated by: 



rime 



(2) 



where Qjet is the jet power, while t) is the jet efhciency and 
we make the assumption that it is only a function of spin. 

Our definition o f the jet ef ficiency is the s ame as that 
of iMerloni fc Hein j (|2008l ) and IShankar et al.l ifioOaV The 
main conceptual difference is that we assign it to black hole 
spin and use the theoretical results from the l iterat u re. W e 
note that our method is different to that of iDalvl l|2009l ). 
since that work assumes the jet power to depend on black 
hole spin and mass, rather than spin and accretion rate. 

It is beyond the scope of this work to discuss in detail 
the physics of jet formation around spinning black holes. We 
use recent jet efhciencies from the literature, and refer the 
reader to that work for the details and assumptions (we list 
these in the following paragraph and in Table [1]). However, 
we do outline some of the key aspects and limitations of 
these efficiencies and how they can affect the results of our 
work. 

The models we will be using are from 
iTchekhovskov et all lj2010h hereafter TIO (we will actually 
use tw o separ ate results: TlOA and TlOB), Benson fc Babul 
(2009', BB09). lNemm en et al' (2007|, N07). lHawlev fc Krolik 
(2006, HK06) and McKinncy (2()0i, M05). In figures and 
tables they will be listed in reverse chronological order. 
Their dependence on the spin a is shown in Figure [2] 

Assuming an efficiency which is solely a function of spin 
is undoubtedly a simplification which sweeps many issues 
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Figure 2. Jet efficiencies ?j as a function of black liole spin, a. The different efficiencies fr om the literature are shown in reverse 
chronological order, so most re cent first. The ac r onym s are: TI P for iTchekh ovskov et al. I ||2010| '). TlOA is their mo d el H/ R=0.5, TlOB 
the H/R= 1 mo del. BB09 for lBenson fc Babull ||2009|') . N07 for lNemmen et al.l l l2007t) , HK06 for lHawIev fc Kroliki l l200d ) and M05 for 
iMcKinnevI ||2005| '). The curve depicting results of HK06 is the analytic approximation found by these authors (their Equation 10), while 
for M05 we have interpolated between the discrete values given by their simulations. 



under the carpet. One such issue is the possible dependence 
of accretion rate with spin, which we assume is small enough 
to have no significant effect (note that the efficiencies from 
simulations, M05 and HK06, will have this effect included). 
Additionally, variations in the magnetic field geometry can 
lead t o similar accretion fiow s producing very different jets 
(M05 iBeckwith et~alll2008al ) , another effect that we are not 
able to take into account here, other than by considering a 
variety of different jet efficiencies. 



Another effect that has been neglected is the effect the 
thickness of the accretion disc will have on the jet power, 
which has been discussed in detail in the li t erature (e.g . 
iGhosh fc Abramowicj[T997l : iLivio et al.|[l999l : iMeieJ l200ll . 
TIG). 

On the one hand, it has been argued that thicker discs 
can sustain stronger poloi dal magnetic fields and hence lead 
to more powerful jets (e.g. lLivio et al.lll999l : lMeiejl200ll ). On 
the other hand M05 and TIO argue that thick discs cover 
a large fraction of the area of the black hole, so that the 
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Shortname 


Reference 


h/r 


TlOA 


Tchckhovskov et al. (2010) 


0.17" 


TlOB 


Tchckhovskov et al. (2010) 


0.33" 


BB09 


Benson & Babul (2009) 


0.35-0.45 


N07 


Nemmen et al. ( 2007) 


0.9'' 


HK06 


Hawlev & Krolik (20061 


0.18-0.4'= 


M05 


McKinnev (20051 


0.2 



Table 1. Summary of the jet models used here. The first and 
second columns state the shortname we will use to refer to the 
efficiency and the bibliographic reference, respectively. The last 
column states the range of accretion flow aspect ratios in the mod- 
els or simulations. "The jet efficiencies for TIO were estimated by 
A. Tchekhovskoy to be r]= P'yv'^ /[2cn{h/r)], where P is the jet 
power shown in Figure 6 of TIO (A. Tchekhovskoy, 2010, priv. 
comm.). This estimate is obtained by assuming that at a radius of 
GM/c^, the magnetic pressure equals the ram pressure in the ac- 
cretion disc, and assuming the radial infall velocity of the accreted 
material in the jet-producing region to be ■yv^oi—c. The thickness 
quoted by TIO includes the corona and jet and corresponds to 3x 
the thickness of the accretion disc. Hence the H/R=0.5 model of 
TIO, labelled here as TlOA, is assigned a scale height h/r=0.l7, 
while the H/R=1.0 model of TIG is assigned h/r=0.33 and la- 
belled TlOB. ''Assuming a viscosity parameter a= 0.3. '^These as- 
pect ratios are taken from lDe Villiers et al. I 1I2OO3I') . 



energy outflowing from this area will go to the accretion 
disc and not into the jet. For this reason, M05 makes a 
diflterence between "total" and "jet" powers. TIO assume the 
magnetic flux through the black hole is constant, not higher 
with larger disc thickness, so that the occultation eflFect leads 
to lower jet powers for thicker discs (TlOB) compared to 
thinner ones (TlOA). The fact that some of the energy can 
go into the disc rather than into powering the jet is also 
mentioned in HK06, although it is not taken into account 
into the efficiency we use (their Equation 10). 

Punsly (2007, 2010) discusses a component of the jet 
driven by the disc in the ergosphere. This component is only 
present in the BB09 and HK06 efficiencies. For HK06, how- 
ever, we are using their analytic approximation which does 
not always reproduce exactly the results from their simula- 
tions. In particular their simulations have maximum spins 
of a=0.99, while we will be using the approximation up to 
a=0.998, a spin for which no simulations have been carried 
out. 

In order to calculate efficiencies, M05 and HK06 used 
the accretion rates from their simulations. TIO present di- 
mensionless jet powers only, and to convert these to effi- 
ciencies one must assume an accretion rate. This accretion 
rate was estimated by assuming that at a radius of GM/(? 
the magnetic pressure in the approximately equals the ram 
pressure of the accretion disc, and assuming the radial in- 
fall velocity to be ~c in this region (A. Tchekhovskoy , 
2010, priv. comm., see Table [T]). The M05 simulations have 
a magnetic pressure that is typically (but not always) lower 
than the ram pressure. For this reason the TlOA efficiency 
is higher than M05, despite having a similar disk thickness, 
/i/r~0.2. The TlOB efficiency, on the other hand, is very 
similar to the M05 one as can be seen in Figure (2] 

Neither M05, HK06 or TIO include cooling in the accre- 
tion disc, while N07 and BB09 assumed an ADAF accretion 



flow to estimate the accretion rate, with negligible cooling. 
Hence, all of these jet efficiencies are estimated for accretion 
flows with no signiflcant radiative cooling. One concern is 
therefore whether the jet powers from these models will be 
appropriate for AGN with high accretion rates. 

Here we justify our reasons for using Equation [21 in 
order to minimise the effect of these uncertainties. 

1. Keeping the jet power an explicit function of both jet 
efficiency and of accretion rate. This will partly account for 
the fact that different spins will lead to different accretion 
rates, by keeping the total jet power explicitly dependent on 
both. It also attempts to reflect the fact that the magnetic 
fields powering the jets origin ate in the accretion di sc and 
are not independent of it (e.g. iDe Villiers et al]|2005l . M05, 
HK06). Given that the accretion rate depends on the black 
hole mass via an Eddington ratio mocAm,, there is a natu- 
ral scaling where black holes with higher masses can achieve 
higher jet powers. However, the jet power is not only depen- 
dent on m. , it also depends on the amount of accretion the 
black hole is undergoing. 

2. Using an efficiency derived for accretion flows with 
negligible cooling. While this appears as a major shortcom- 
ing for black holes with high Eddington rates, we believe it is 
still a useful approximation. The higher accretion rates are 
already considered in the m term. The worry is then if the jet 
efficiency is appropriate for an accretion flow with a high A. 
As mentioned earlier, the aspect ratio is believed to be an im- 
portant factor in determining the jet power. The jet efficien- 
cies used here are derived from accretion flows with typical 
aspect ratios 0.17^/i/r^0.9, therefore covering a wide range 
of disc thickness. In this paper we will concern ourselves 
mainly with two types of AGN: those at very low Edding- 
ton ratios (ADAFs) and those at very high ones (quasars). 
Equation[2]is appropriate for ADAFs, but an important fact 
is that the accretion flows in quasars are radiation-pressure 
dominated, not gas-pressure dom inated, and will have typ - 
ical aspect ratios hlr ^X ('NT73. [Abramowicz et al.lll988l ). 
For quasars, A~0.25 jMcLure fc Dunlodl2004l )" so that the 
aspect ratio /i/r~0.25. Hence the thickness of the two types 
of flow are similar. This is the reason we make the assump- 
tion that the jet efficiencies r] shown in Figure[2]can be used 
for both quasars and ADAFs. 

It is immediately clear that r) can differ significantly 
from model to model, which will affect quantitatively our 
results. However, given that in all of these cases 77 increases 
with increasing spin a, our results will be qualitatively the 
same for different assumed jet efficiencies. 

We stress that all the jet efficiencies have assumptions 
and limitations. We believe that showing the results for six 
different jet efficiencies from the literature is the most infor- 
mative way of presenting our results. Now that our assump- 
tions have been clearly stated, we proceed to compare the 
results to observations. 



3 THE RADIO LOUDNESS OF QUASARS 

A necessary test of our parametrisation is whether it can 
produce the range in radio loudness observed amongst 
quasars. 

Quasars are SMBHs with black holes of m,^ 
10* M0 accreting close to the Eddington limit (e.g. 
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iMcLure fc Dunlod 120041 '). For a given optical luminosity, 
a proxy for bolometric luminosity in unobscured objects, 
quasars are found to have a vast range of radio luminosi- 
ties. Most of the quasars are found at the 'radio-quiet' end, 
where the ratio of radio to optical luminosity density is ^10. 
However, a small fraction show radio luminosities up to 1000 
times larger, and it has been suggested that spin is the hid- 
den variable dorninating this scatter (e.g . IScheueij 1 19921 : 
IWilson fc Colbert! ll995l : ISikora et al.l l2007l i . 



3.1 Estimating observable quantities 

We use Equations [T] and [2] to predict the observed properties 
associated with quas ars, and test the s e aga inst the sample 
of quasars studied bv lCirasuolo et al] (|2003l '). The accretion 
rates for our predictons are given by: 



(3) 



where Le is the Eddington-rate luminosity, 1.3 x 10''^ W per 
solar mass, m. is the black hole mass, and the bolometric 
luminosity Lboi is given by Lboi = L^Xm ,. Typical values for 
quasars are m, > 10* Mq, A ~ 0.3 (e.g. iMcLure fc Dunlod 
|2004 ). so that Lboi~ 4 x 10^* W and m ~0.7 Mq yr'^ 

To convert from bolometric luminosity, Lboi, to B-band 
monochromatic luminosity, L^b, we use: 



(4) 



where, in the B band, z/b is the frequency and C^b is the 
bolometric correction that converts from the monochromatic 
luminosity to bolometri c luminosity. Here w e use the bolo- 
metric corrections from I Hopkins et al.l (|2007l ). 

To convert Q. et to L^isi we use the conve rsion given by 
iMiller et ai] (| 19931 ') for radio-quiet quasars and lWillott et all 
( 1999) for radio-loud quasars: 



Kw ) 



3 X lO^V^ ( 



1028 WHz~isr-i 



(5) 



The term / represents the combination of several uncer- 
tainty terms when estimating Qjet from L^isi, which in- 
cludes for example the filling factor in the lobes, departure 
from minimum energy, and the fraction of energy in non- 
radiating particles. 

The plausible range of values for / has been estimated 
to be 1 s; / 20 by Willott et al. (1999, for details see Sec- 
tion 4 of their work). An independent method of estimating 
jet power from the radio luminosity has been developed by 
estimating the e nergy stored in X-ray cavities in the intra - 
cluster medium jBirzan et al. I l2004l : ICavagnolo et al.|[201ol ). 



ICavagnolo et al. I l|2010l ) find that /~20 is required for both 
methods to agree. 

For the remainder of this paper, we will therefore use a 
value of f~20, which means the jet powers estimated using 
the re lations of "Willott et al. (1999) and ICavagnolo et al.l 
l|2010h agree with in their unce r tainti es. However, we note 
that t he studies of lBirzan et al.l (|2004l l and lCavagnolo et al.l 
l|20ldl consisted mostly of sources of moderate radio lumi- 
nosity. Sources with higher radio luminosities might have 
typically lower values of /, close to 1, as we discuss below. 

The value of /=20 assumes a large contribution from 
non-radiating particles (~200 times more energy in non- 
radiating particles than in radiating ones, see Section 4 of 



Willott et al. 1999). A large fraction of non-radiating parti- 
cles can be inferred observationally if the pressure provided 
by radiating particles inside the lobes (derived from min- 
imum energy arguments) is much lower than the external 
pressure acting on the lobes due to the intra-cluster gas. 
In this case the extra pressure must be provided by par- 
ticles that are not emitting synchrotron radiation. There 
exist cases where the pressure from radiating particles is 
comparable to the external pressure, so that the fraction 
of non-radiating parti cles is small. In tho se cases /~1 is 
more appropriate (see lCroston et al]|2008l . and Figure 2 of 
Cavagnolo et al. 2010). From studies of the cavities in the 
intra-cluster medium, it seems possible that the sources that 
are found to have low jet powers for a given radio lumi- 
nosity have a lower content of non-radiating particles, due 
to a lower ent rainment of external material by the jets. 
ICroston et al.l 2008: Cavagnolo et al. 2010). 

Fewer detailed studies exist of the more powerful radio 
sources, since these are rare in the nearby Universe. How- 
ever, some studies suggest the fraction of non-radia ting par- 
ticles is small in these cases (|Croston et al.ll2004l ). This is 
possibly due to th e fact t hat the most powerful radio sources, 
iFanaroff fc RilevI (|l974l ) class II objects, have strong bow 
shocks that push back the intergalactic medium so that en- 
trainment is minimal. In these obj ects, a value of f^l might 
be more appropriate in this case (|Croston et al.l l2004'). The 
arrows shown in Figure [3] show how the curves would move 
if / were changed from /=20 to /=1. 

Source-to-source variations are likely important, as re- 
flected by the scatter between jet power and radio lumi- 
nosity (Birzan ct al. 2004; Cavagnolo et al. 2010). However, 
there might be systematic variations with luminosity, with 
the high-radio luminosity sources having a lower fraction of 
non-radiating particles, and hence requiring a lower value 
of /. In such cases, assuming /=20 will systematically un- 
derestimate the radio luminosity that can be produced for a 
given jet power. For a given jet power, L^isi oc / t , mean- 
ing that changing / from 20 to 1 would increase the radio 
luminosity by practically 2 dex (see the arrow on Figure [3]). 

Since the comparison data is at 1.4 GHz, we convert 
from 151 MHz assuming a single power law L^, oc i/"" as- 
suming a spectral index of q = 0.75. 



3.2 Comparison of the models to the data 

Figure [3] shows the loci of radio luminosity density versus op- 
tical luminosity density from our assumptions in Section [2] 
and Equations O |4] and [5] Each panel assumes the jet effi- 
ciencies of Figure [5J in the same order. We have plotted two 
tracks for each jet efficiency: one for a=0.1 (dotted line) and 
another one for a=0.998 (dashed line). The cutoff maximum 
spin is somewhat arbitrary, but our choice is motiv a ted b y 
the theoretical maximum value predicted bv lThornd (jl974h . 

The black points are the data for the quasars studied by 
Cirasuolo et al. (2003), who matched two samples of opti- 
cally selected quasars (the 2dF QSO red shift survey and the 
large bright quasar survev lCroom et"~al. , 2001: Hewett et al.l 
Il995l ) to the FIRST radio survev (|Becker et al.lll995l '). The 
sensitivity of the FIRST survey is low enough that several 
decades of radio loudness can be explored, yet the sample 
is probably incomplete below radio luminosity densities of 
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Figure 3. Radio vs optical luminosity densities for a sample of z ^2.1 quasars. These were selected by matching an optically- 

selected sample of quasars to a survey at 1.4 GHz (Cirasuolo et al. 2003). Only sources with large enough radio luminosities to make the 
radio flux density cut make it into the sample, so that it is not complete below ~ 10^^ W Hz~^ sr~^. The accretion rates are estimated 
from Equation [3] for the range 3 X 10^ ^Am,^ 3 X 10^ M©. The dotted and dashed lines show the resulting tracks for spins a=0.0 and 
0.998 respectively, using the jet efficiencies rj plotted in Figure[2] The cross shows the 95% confidence interval for the assumed bolometric 
corrections (horizontal hne) and the effect of varying the radio spectral indices in the range 0.25 ^ « ^ 1.25 (vertical hne). The arrow 
shows the effect on the tracks of changing / from 20 to 1. The lower dash-dotted horizontal line repre sents the radio luminosity from a 
star burst forming massive stars at a rate of 1000 Mq yr~^ (for stars more massive than Ss5 Mq only, ICondon|[T993) . This corresponds 
to a total star-fo rmation rate of ^^4000 Mq yr~^ assuming a Salpeter initial mass function. The upper dash-dotted line shows the 
FRl/FRll divide jPanaroff fc Rilevlll974l'l . 



^'1^1.4 ~ 10^* W Hz ^ sr '^ (as noted by Cirasuolo et al. 
2003). 

The data show three main features: a high-radio lumi- 
nosity envelope, a low-radio luminosity envelope, and as a 
result of the difference between these two, a scatter of ~3 



dex. The high-luminosity envelope is not a selection effect 
but the lower envelope might well be due to incompleteness 
below ~ 10^* W Hz-i sr'^ 

Given a maximum achievable jet efficiency, 77(a~l), 
Equations [2] and [1] do predict an upper envelope showing a 
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correlation between radio and optical luminosity, Qjot oc rh. 
The exact location of the envelope, however, depends on 
the model assumed for rj, and different models reproduce 
the data in Figure [3] with varying degrees of success. 

The BB09 model and H[K06 produce upper envelopes 
in reasonable agreement with the data (although we stress 
the caveats mentioned in Section [2]2]). The agreement is not 
perfect and some sources appear ~1 dex higher than the 
curves. However, given the uncertainties in converting from 
bolometric to monochromatic luminosity, from jet power to 
radio luminosity density and the uncertainties in the actual 
radio spectral indices, we do not consider this discrepancy a 
worry. 

The TlOA and N07, efficiencies produce upper en- 
velopes that are significantly lower. This is because the max- 
imum efficiency these two models reach is r;~0.2-0.3, while 
the most powerful radio-lo ud quasars and radio galaxies 
require jet efficiencies ^1 jPunslvl l2007l : iMcNamara et af] 
l201ll : lFernandes et al.ll201oh . Therefore, any model or simu- 
lation that does not reach such high valu es will not be able to 
reproduce the upper envelope (see also |Punslvll201ol ). This 
is the case for the TlOB and M05 models, which produce 
even lower tracks. 

However, we remind the reader that in Section [3.11 we 
have assumed a value of f=2Q when converting from jet 
power to radio luminosity density. This implies a large com- 
ponent of non-radiating particles. The arrows shown in Fig- 
ure [3] show how the curves would move if / were changed 
from /=20 to /=1. As discussed in Section [3. II it is possible 
that the most powerful radio sources systematically require 
lower values of /. For f—1, the high-spin tracks for TlOA, 
N07 and HK06 will easily produce high enough radio lumi- 
nosities, the BB09 will predict far too much radio luminosity 
for a~l, and the TlOB and M05 high spin tracks will be high 
enough to match the upper envelope. 

The upper enve l ope c an also explain the result of 
I Rawlings fc Saunders! l|l99ll ). who found a linear correla- 
tion between jet power and narrow-line luminosity (Lnl, 
an isotropic proxy for Lboi) for radio-loud AGN (their work 
assumed /~2 to convert from radio luminosity to jet power). 
If we assume Lboi ~ 500 — 1000 x Lnl, our the curves with 
maximum effi ciencies r;~0. 1-1.0 can also explain the data of 
iRawlings fc Sa unders ( 1991) as an approximate upper enve- 
lope. ISerieant et al.l (|1998) found an upper envelope for the 
radio luminosity density versus optical magnitude of steep- 
spectrum quasars, and argued it was not due to any se- 
lection effects or relativistic beaming. A similar envelope is 
also found between infrared luminosity density (another ap- 
proximately isotropic estimator of Lboi) and low- frequency 
radio luminosity density for radio-selected AGNs, an enve- 
lope that can again be reasonably reproduced by our curves 
l|Fernandes et al.ll2010l l. 

The minimum jet efficiency, ri{a=0), from HK06, N07 
and BB09 predict lower envelopes, which appear in good 
agreement with the data. However, as Cirasuolo et al. (2003) 
point out, the sample of quasars might be incomplete below 



The fact that in the the M05, TlOA and TlOB mod- 
els as a— >0, means that they do not have a lower en- 
velope, but this is not in disagreement with current data, 
as explained above. The existence of a possible lower enve- 
lope requires further studies using deep radio observations 
of powerful quasars. 

Discriminating whether this lower envelope is a real 
physical effect or just a selection effect, is further com- 
plicated by the fact that very powerful starbursts, with 
massive star- formation rates ~1000 M0 yr~^ (for stars 
^5 Mq, which would correspond to a total star-formation 
rate of ~4000 Mq yr~^ for a Salpeter initial mass func- 
tion), can show radio luminosity densities as high as ^ 
3 X 10^^ W Hz"^ sr"^ without the presence of an AGN 
(this is s hown by the bottom horizontal line in Figure |31 
following ICondon|[l993) . 

3.3 Implications for AGN 

Also overplotted on Figure [3] is the divide between Fanaroff- 
Riley class I and class H objects (jFanaroff fc RilevI Il974l . 
hereafter FRI and FRH sources) . FRI and FRH sources have 
powerful extended jets, and the bulk of the quasars in the 
figure lie around or above the FRI/FRH divide. This sug- 
gests that extended jets are common around quasars, and 
indeed the common assumption that quasars were never 
associated with FRI-type jets has been proven false (see 
iBlundell fc Rai lings 200l|; Ikevwood et~ani2007l l . 



Miller et al.l (| 19931 ') also found that many of the bright- 



L^iA ~ 10^"* W Hz ^ sr ^, so that the agreement between 5 x 10^* W Hz ^ sr 



est optically-selected quasars at 0.3 ^ z ^ 0.5 were FRH 
sources, and the high jet powers inferred bv lCao fc Rawlings! 
(j2004 ) for powerful nearby FRIs (^10^*^ W) seem difficult 
to reach without high accretion rates. 

Turning this around, even with a high spin, no AGN 
with a low radiative luminosity is expected to show FRH lu- 
minosities (and presumably by extension no FRII jet struc- 
tures either). For example, the FRI/FRII division is at a 
151 MHz radio luminosity of ~ 2 x 10^^ W Hz'^ sr'S 
which corresponds to Qjct~ 1-3 x 10^** W (assuming /=20). 
Even with jet efficiencies as high high as r;~0.5, this still re- 
quires an accretion rate of 0.05 M© yr~^, which can only be 
achieved by SMBHs with Am.> 3 x lO'' Mp . This is the typ- 
ical v alue for optically-selected quasars ijMcLure fc Dunlod 
|2004 ). and corresponds to moderate to high Eddington rates 
for black hole masses ~10^ — 10^ Mq, and only corre- 
sponds to modest Eddington rates for the rare black holes 
with masses ^3 x 10^ Mq. This explains why FRII sources 
are mostly associated with HEG spectra, and only some- 
times with LEGs spectra, while FRIs whi ch can be associ- 
ated with both HEG or LEG spe c tra (e.g. Ledlow fc OwerJ 
! 19961 : !Blundell fc Rawlings! !200ll : !Hardcastle et al.! !2006! ). 
Section [5] describes this quantitatively. 

An upper limit on the observed radio luminosity den- 
sity, 1/1,151, can be estimated by assuming 7)max~l, emax~0.3 
as well as Amax ~ 1 and m.,max~ 3 x 10^ Mq. This yields 
maximum values of Qjct,max ~ W, and Li^isi,, 



i jct ,max 

or L^i.z 



1 X 10^^ W Hz"^ sr" 



the lower envelope and the data might be fortuitous. Again, 
the uncertainties when converting iboi and Qjot into L^b 
and I/I/151 are large (see Figure so that quasars with 
slightly lower radio-to-optical luminosity density ratios are 
not strictly in disagreement with our curves. 



No radio galaxy observed so far has a radio luminosity sig- 
nificantly above th is value (e.g. !Rawlings fc SaundersI 199l! : 
!Willott et al]!l999! ). 

We now proceed to test whether under reasonable as- 
sumptions the spin paradigm can explain the different con- 
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tributions to the local radio LF from AGN with high and 
low accretion rates (HEGs and LEGs, respectively). 



4 MODELLING THE LOCAL RADIO 
LUMINOSITY FUNCTION 

We proceed to fit the local radio LF of LEGs and HEGs 
to find the distribution of spins that best describes each 
subpopulation. The best-fitting parametrisation allows us 
to interpret the physical conditions that lead to the spin 
distribution, it can be used to derive an estimate of the spin 
history of SMBHs. 

Above ^ 10^^ W Hz^^ sr-\ all the sources 

contributing to the radio LF have m.^ 10* MQ(e.g. 
iMcLure et aLll2004 ISmolcic et alll2009l ). We will therefore 
limit our discussion to SMBHs with m.^ 10* M©. The ra- 
dio luminosity is determined from the jet power Qjct (see 
below), which from Equation[2]depends on the jet efficiency, 
T) as well as the accretion rate, m. 

Throughout this section we follow the convention that 
the sign <^(L) refers to a luminosity function while the sign 
P(x) refers to a normalised probability distribution function 
(PDF): 



P{x)dx = 1 



(6) 



Luminosity functions are PDFs that are not normalised to 
unity, but rather to a space density. Hence we will manipu- 
late LFs in an analogous way to PDFs (note that in practice 
all the integrals are done using base- 10 logarithms). 

Assuming that the jet efficiency is essentially a mono- 
tonic function of spin, the PDF for jet efficiencies is related 
to the PDF for spins by: 



P{rj)dri = P{a)da 



(7) 



The space density of sources with a given radio lumi- 
nosity density is related to the space density of sources 
with a given jet power Q, via: 



<f>{L^)AL, = <f>{Q)AQ. 



(8) 



From Equation [21 the distribution of jet powers will 
depend both on the distribution of accretion rates and spins. 
The PDF of a given jet power Q can be thought of as the 
probability of Q and an accretion rate m and an efficiency 
77, marginalised over the 'nuisance' parameters fn and rj. 



dQ, 



jet 



PiQ. 



dQ. 



jot 



dL„ 



di, 

P{Qjct,m,rj)dmdr] 



(9) 
(10) 



Using Bayes' theorem, we can write this in terms of the 
PDF of Q given m and rj and the PDFs of rh and rj. 



P{Qjct,m,r]) = P{Qjct\m,ri)P{r],m) 



(11) 



For a given accretion and jet efficiency, only one value 
of the jet power is possible, so that the probability of Qjct 
given rh and 77 is a delta-function: 



-P(Qjctlm, 77) = S{Q jet — rime ), 



(12) 



and hence 



P{Qiat\m,ri)P{r),rn)drhdri (13) 
5(Qjct — 77mc'^)P(r;, m)drhd77 

= j P{rj,rh' = Qjct - r]c )dr) 

For a particular class of objects, we make the approxi- 
mation that the spin is independent from the instantaneous 
accretion rate, so that: 



P{r),rh') = P{ij)P{rh'). 



(14) 



The space density of sources with a given radio lumi- 
nosity is therefore given by: 



0(L,) 



^ / P(,7)0(m )dr; 



dQjct 
dL,. 



P[a)(f>(rn')da 



(15) 



We model the LEGs and HEGs separately, and we begin 
with the HEGs. 



4.1 High-accretion rate radio galaxies 

The HEGs are modelled by SMBHs with high accretion 
rates, which we refer to as QSOs for simplicity. The space 
density of SMBHs with a given accretion rate can be esti- 
mated from the X-ray LF, since the X-ray luminosity can 
be used to estimate the bolometric luminosity via 



(16) 



where Cx is the bolometric correction at X-ray energies. 
Hence, the accretion rate m can be estimated as: 



CxLy 



(17) 



The space density of HEGs with a given radio luminos- 
ity is therefore: 



<^QSo(ii^) = 



dQ. 



jot 



dL„ 



J PQSoi&)<t>{L'x)^da 



(18) 



To avoid including SMBHs with m,< 10* M0, we im- 
pose a minimum bolometric luminosity of Lboi ^ 3 x 10'^* W. 
This can only be achieved by SMBHs with m.> 10* Mq at 
25% of t he Eddington luminosity (the characteristic value 



Ao/o 01 t ne Ji^ddmgton luminosity 
found bv lMcLure fc DunlodlioM ). 



Assuming typical bolometric corrections, this corre- 
sponds to an X-ray luminosity of 9 x 10''^ W. We have 
also limited the upper X-ray luminosity to ^ 1 x 10'^^ W 
since no sources brigh ter than this have been observed (e.g. 
[Silverman et al.ll2008l ) 

The X-ray LF of [Silverman et all (|2008l ) is selected in 
the hard band (2-8 keV), meaning it is sensitive to unab- 
sorbed and moderately absorbed sources (with column den- 
sities A^H^IO^'^ m"^ and lO^** m"^ sJA^hsS 10^* m"^, respec- 
tively). However, sources with column densities ^1.5x10^* 
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are optically-thick to Compton scattering ('Compton- 
thick') and will be missed by the X-ray LF. 

In the X-ray luminosity range of interest, 9 x 10^® W 
<Lx^ 1 X 10^^ W, approximately half of the hard X-ray 
sources are unabsorbed and the other half are absorbed but 
Compton thick. The fraction of Compton-thick AGN is es- 
timated to be similar to the fraction of Compton-thin ab- 
sorbed AGN, meaning the space density of quasars given 
by the hard X-ray LF has to be corrected by a factor ~1.5 
(see e -K- lRisaliti et al.|[T999l : lGim et al.ll2007l : fGoulding et al] 
l20ld '). Although there migh t be a slight evolution in the frac- 
tion of absorbed AGN (e.g. iTreister fc Urrvll2006l : iHasinged 
l2008h . thi s correction is likely still cor r ect at high red- 
shift fe.g. [Marti'nez-Sansigre et all l2007l : iGilli et all l2007l : 
' 20091 ). We therefore multiply the space density 



iFiore et al 



of the X-ray LF by a factor of 1.5, to estimate the space 
density of all QSOs (whether hidden or not). 

4.2 Low-accretion rate radio galaxies 

We model the LEG component of radio LF as SMBHs with 
low accretion rates, which we call ADAFs. This requires 
some extra assumptions, since no luminosity function is 
available to estimate their accretion rates. However, the local 
mass function of supermassive black holes, hereafter BHMF, 
provides the distribution of masses, which can be related to 
the distribution of accretion rates via a distribution of Ed- 
dington rates: 



P(m') = j P(m ,m,)dm, 
= J P{m' \m,)P{m,)dm, 



(19) 



We have expressed the conditional probability of an ac- 
cretion rate given a black hole mass, P{rh'\m,), as a proba- 
bility of an Eddington rate -P(A). 

The radio LF due to ADAFs can therefore be modelled 

as: 



4>AUAFiL^) — 

dQjet f p^^^p(^) j' p(A)0(m.)^dm.da 



dL, 



(20) 



Given our ignorance about the distribution of Edding- 
ton ratios, we choose a prior PDF which is flat in logarithmic 
space, known as a Jeffreys prior: 



P{\) oc 



(A„ 



Amin) 



(21) 



We assume the ADAFs are described by SMBHs accret- 
ing at Eddington ratios in the range 10"** s;A< 10-^ hence 
A„,ax=10~2 and An,in=10-^ 



4.3 Model selection for the spin distribution 

We now proceed to use the data from the observed local ra- 
dio LF of LEGs and HEGs to constrain the spin distribution. 
Three models are considered. 



Model A consists of a power law PDF for the spin: 

Pa (a) = NAa-" (22) 

The power law index p is the only free parameter in this 
model, and it is given a flat prior PDF in the range —4 ^ 
p ^ 4. Na is simply a normalisation constant. 

The next model, B, is a single gaussian distribution: 



1 ( f»£ 
-2( c 



(23) 



We fix (7a=0.05 and Nb is the normalisation term, so this 
model has also only one free parameter, fiB, to which we 
assign a fiat prior in the range ^ /is ^ 1. 

The value of CTa=0.05 is assumed since it can provide a 
reasonable approximation to the width of the spin distribu- 
tions expected to arise from accretion onto SMBHs as well 
as from major mergers (see Section [S]). Note that if /is is 
close to one edge, the PDFs will be truncated, so we calcu- 
late the normalisation term numerically, rather than using 
the analytic value of A'^b = l/aaV^iv. 

Model C consists of two gaussians distribution with 
variable relative height: 



Pc(a) = A^c 



Hcie 



+ e 



(24) 



There are three free parameters: the two mean spins fici, 
^C2 as well as the relative height of the first gaussian Hci- 
The standard deviation of both gaussians is kept constant 
at the value (Ta=0.05, and A'^c is a normalisation term. All 
three parameters, /ici and ^C2 are assigned fiat priors in 
the range 0-1, while due to our ignorance of scale, Hci is 
assigned a fiat prior in logio, in the range -3 to 0. 

For each model, we calculate the 'likelihood' 
P(data|{2;i, ...Sn}, model) as a function of the n-parameters, 
{2:1, ...Xn}: 



PADAF(data|{a;i, ...In}, model) = 
Padaf (datai | { xi , . . . Xn } , model) = 



,(. 



iLEG -'^'ADAF ) 
''■LEG / 



'^<*LEG v27r 



(25) 



and 



PQSo(dataj{2;i, ...Xn}, model) 
Pqso (datai I { Ki, ...Xn}, model) 

i 

1 '1. 



'^</>HEGV27r 



(26) 



where ^iLEG and ^iheg are the data points from the ob- 
served radio LFs for LEGs and HEGs respectively, while 
0ADAF and 0QSO are the modelled radio LFs of high- and 
low-accretion rate SMBHs. These are independent, so that 
the combined likelihood is given by: 



Ptotai(dataj{2:i, ...Xn}, model) = 
PQSo(data|{xi, ...Xn}, model) Padaf (dataj{ }, model) (27) 
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Models 


logio[Odds Ratio] 


TlOA 


B / A 


257 


TlOA 


C / B 


100 


TlOA 


C / A 


357 


TlOB 


B / A 


>247 


TlOB 


C / B 


34 


TlOB 


C / A 


>281 


BB09 


B / A 


>177 


BB09 


C / B 


31 


BB09 


C / A 


>209 


N07 


B / A 


>261 


N07 


C / B 


51 


N07 


C / A 


>311 


HK06 


B / A 


>161 


HK06 


C / B 


3 


HK06 


C / A 


>163 


M05 


B / A 


>264 


M05 


C / B 


54 


M05 


C / A 


>319 



Table 2. Summary of results of the model selection. Section 14.31 
describes the models and the model selection method. The first 
column indicates the two models being compared, with the most 
likely model on the left. The second column shows the base-10 
logarithm of the odds ratio (given by Equation (30] see Section [4.3l 
for more details). The right-censored results marked with a > 
occur due to the fact that the evidence of model A is <10^^^*, 
too small for the program. 

The integral of the likelihood over the parameter space 
{xi, ...Xn} yields the 'evidence', P(model|data): 

P(data| model) = 
/ P({xi, ...Xn})Ptotai(model|{ },data)d"x, (28) 

and Bayes' theorem allows us to infer the probability of the 
model given the data: 

P(modelldata) = J^(data|model)P(model) 
^ ' ' P(data) ^ ' 

Where P(model) is the prior for a model. When comparing 
two models, the denominator can be eliminated by taking 
the 'odds ratio', P(modell|data)/P(model2|data). 

P(modell|data) _ P(data|modell)P(modell) 
P(model2ldata) ^ P(datalmodel2)P(model2) ^ ' 

Hence, using the odds ratio, we can compare models 
two at a time, and keep the model which fares best agains 
the others. In Equation 1281 the integral over the parameter 
space, and the priors P({xi, ...Xn}) penalise the use of extra 
parameters as well as the exploration of parameter space 
with a low likelihood. Hence, models with different number 
of parameters can be compared in a fair way. 

Finally, we use the different jet efficiencies 77(a) de- 
scribed in Section [5] and shown in Figure [51 Using these 
jet efficiencies and our parametrisation we are able to fit the 
local radio LF, and Section [5] presents the results. 



V 


Population 


MCI 


P-C2 


Hci 


TlOA 
TlOA 


ADAF-LEGs 
QSO-HEGs 


0.99 
0.70 


0.29 
0.01 


0.316 
0.079 


TlOB 
TlOB 


ADAF-LEGs 
QSO-HEGs 


0.99 
0.99 


0.99 
0.00 


0.398 


BB09 
BB09 


ADAF-LEGs 
QSO-HEGs 


0.89 
0.00 


0.08 
0.00 


0.631 


N07 
N07 


ADAF-LEGs 
QSO-HEGs 


0.99 
0.99 


0.16 
0.00 


0.501 
0.040 


HK06 
HK06 


ADAF-LEGs 
QSO-HEGs 


0.99 
0.00 


0.06 
0.00 


0.794 


M05 
M05 


ADAF-LEGs 
QSO-HEGs 


0.09 
0.99 


0.99 
0.09 


0.100 
0.398 



Table 3. Summary of best-fitting parameters for model C, chosen 
ahead of the other two models. The quasars (QSOs) have been 
modelled using the X-ray LF as described in Section l4.1l and the 
parameters are those that best fit the observed radio LF of HEGs. 
The ADAFs are modelled using the local BHMF, as described in 
Section |4. 21 and have been compared to the observed radio LF of 
LEGs. 

5 THE DISTRIBUTION OF SPINS FROM THE 
LOCAL RADIO LF 

5.1 Comparison to the z=0 radio LF 

To constrain the best-fitting model and parameters, we use 
the data kindly provided by P. Best prior to publication 
(Best et al. 201 1, in prep., from the original sample of 
iBest et aljlioosi ). who have separated the local radio LF 
into HEGs and LEGs. These are plotted as empty triangles 
and filled circles, respectively, in Figure [l] 

Using Equation 1301 and assuming equal priors for all 
models, P(modelA)=P(modelB)=P(modelC), we select the 
best model, which we find to be model C (see Table [2]). 
We note that model C has been chosen despite the Ockham 
factors that penalise the exploration of extra parameters. 

Once we have chosen a model, the best-fitting parame- 
ters are chosen by maximising the likelihood, which is equiv- 
alent to the traditional minimisation. We find the best- 
fitting parameters for SMBHs with high and low accretion 
rates ('QSOs' and 'ADAFs' respectively) separately, and Ta- 
ble |3] summarises the results. Table |4] summarises the best- 
fitting parameters for the rejected models. 

Figure [4] shows that the best- fitting parameters can ac- 
curately reproduce the local radio LF of LEGs and HEGs 
in the region ^ 3 x 10^^ W Hz"^ sr"^ The region 

1-3x10^'^ W Hz~^ sr~^, is not so well fit, but a few points 
are worth mentioning. Although the fits look better if we 
only fit the region > 3 x 10^^ W Hz"^ sr"\ the best-fitting 
parameters are essentially unchanged. 

Some models for the efficiencies do a remarkable job of 
explaining the radio LF, particularly the LEGs (e.g. BB09, 
HK06) . Models that cannot reach very high values of r] strug- 
gle to fit the LEGs, but can still fit the HEGs reasonably 
well. For example the shape of the QSO model function from 
M05 is almost identical to the HEG LF, although the nor- 
malisation is slightly off. 
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Parameters for model A 





P ADAF 


P QSO 




(a>QSO 


{o)sMBH z=0 


{o)sMBH z=l 


TlOA 


0.5 


1.3 


0.34 


0.05 


0.32 


0.20 


TlOB 


-4.0 


0.7 


0.83 


0.26 


0.42 


0.27 


BB09 


0.3 


4.0 


0.41 


0.00 


0.25 


0.04 


N07 


-0.6 


4.0 


0.61 


0.00 


0.52 


0.16 


HK06 


-3.2 


4.0 


0.81 


0.00 


0.59 


0.12 


M05 


-4.0 


0.5 


0.83 


0.34 


0.71 


0.42 


Parameters for model B 














'n 


Mb ADAF 


MB QSO 


{o)sMBH z=0 


{a)sMBH z=l 






TlOA 


0.48 


0.13 


0.32 


0.16 






TlOB 


0.04 


0.49 


0.49 


0.49 






BB09 


0.85 


0.04 


0.61 


0.15 






N07 


0.79 


0.04 


0.53 


0.12 






HK06 


0.91 


0.04 


0.66 


0.16 






M05 


0.96 


0.73 


0.92 


0.78 







Table 4. Summary of best-fitting parameters for the rejected models: A and B. The model selection and parameter-fitting procedures 
are described in Section 14.21 Models A and B show the same trend of evolution in {a{z)) with z as seen in model C. 



We have assumed a value of /=20, but if we change this 
to f=l (at all radio luminosities) the TlOB and M05 models 
reproduce the LEGs better, and provide a reasonable fit. In 
this case the M05 model requires all the spin to lie around 
a=0.14 for the HEGs, while the LEGs require a bimodal 
distribution (with means of /ici=0.86 and /ie2=0.10). The 
TlOB requires bimodal distributions for both the HEGs 
(0.33 and 0.00 respectively). The TlOA, BB09, N07 and 
HK06, on the other hand, require all the probability den- 
sity to lie around zero spin (both /ici and /ic2=0.00 for 
both HEGs and LEGs). Despite the low spin values, all three 
models overpredict the LEGs and HEGs, and provide a very 
poor fit in the case of f=l. 

As was discussed in Section 13.11 the typical value of / 
might change systematically with luminosity. A possible im- 
provement to our modelling would be to vary / with radio 
luminosity. Currently we do not wish to follow this approach 
since we do not have enough independent constraints, so al- 
lowing / to vary constrained only by the radio LF of HEGs 
and LEGS would lead to 'fine tuning' of parameters rather 
than allowing us to make inferences. The value of f—20 
is justified independently and is appropriate for the radio 
sources in the local Universe (see Section [5TT1) . 

The QSO LF includes a break (due to the break in the 
X-ray LF) which looks similar to the break observed in the 
HEGs, although the normalisation is not exactly in the right 
place. The shape of the break in the X-ray LF is approxi- 
mately preserved due to the fact that the jet efficiencies 
are generally shallow at the low-spin end (e.g. HK06, N07, 
BB09). The jet efficiencies that are steeper at the low-spin 
end have smoother breaks (M05, TIO). 

Similarly, in many models there is a break in the LEG 
LF. This originates in the break of the mass function, but 
in some models with very steep values of i] at high spin the 
break is smoothed out. This is most noticeable in BB09. 
The observed break in the LEG radio LF is not as marked 
as that in the mass function since, due to their different 
accretion rates and jet efficiencies, SMBHs with a range of 
masses contribute to that region of the LEG radio LF. 



We note that removing the Compton-thick correction, 
modelled as a 50% increase in the space density of HEGs 
(Section 14. ip . makes no difference to our results. For all the 
jet efficiencies, the resulting best-fitting parameters for the 
HEG-QSOs are practically identical to those when including 
the Compton-thick sources. The only exception is for BB09: 
when removing the Compton-thick correction the smaller 
gaussian has a mean spin of 0.99 and an amplitude 0.01 
relative to the gaussian with mean spin 0.00 (compare with 
the results in Table [S]). 

Overall, we consider the fits acceptable, given that they 
are based on a few simple physically-motivated prescrip- 
tions. We remind the reader that the uncertainties in con- 
verting jet powers to radio luminosity densities, discussed in 
Section [3.11 applv here equally. 

Most notably, assuming /=20 requires the models to 
achieve ?7max~l in order to explain the LEGs, due to their 
low accretion rates but high radio luminosities. Reducing 
/ to ~1 would mean the low- 77 models would fare much 
better (e.g. M05, TlOB and to a certain extent N07) but 
then the models with high rymin, such as HK06 and BB09, 
would overpredict the radio LF. 

The exact choice of cutoff for the maximum spin also 
makes a difference. For example BB09 argue for lower values 
of maximum spin in ADAFs, and indeed if we set an upper 
cutoff of a~0.95 the BB09 model provides a better fit to the 
LEGs than for o=0.998. 

We have only fitted the data above radio luminosity 
densities of 10^^ W Hz"^ sr"\ but we can see how our best- 
fits fare in the lower luminosity region, where we have not 
fitted to the data. Figure [S] shows what our models predict 
in this region. Indeed, extrapolation of our model radio LFs 
to lower radio luminosity densities does not overpredict sig- 
nificantly the observed radio LFs of HEGs or LEGs. 

The ADAF models provide a very good match to the 
LEGs in the range 10^^-10^^ WHz~^ sr"\ There is no need 
for a further component to explain the LEGs: the prediction 
from our model is therefore that SMBHs with m.^ 10* Mq 
dominate the LEG radio LF down to 10^^ W Hz'^ sr"^ 
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Figure 4. The z=0 radio luminosity function split into iiigh-excitation gas (HEGs, open triangles) and low-excitation gas sources (LEGs, 
black circles, data from P. Best, et al., in prep.). The empty squares indicate the summed radio LP. The acronyms summarise the jet 
efficiencies 77 shown in Pigure[2] The curves show the best-fitting curves from our most appropriate model (model C), chosen using the 
odds ratio as described in Section l4.3l The dashed line shows the radio LP modelled from high-accretion rate SMBHs (labelled QSOs, for 
simplicity), and is compared to the HEGs, while the dotted line shows the radio LP modelled from low-accretion rate objects (labelled 
ADAPs) which is compared to the LEGs. The model has three free parameters which describe the spin distribution. These parameters 
are fixed by comparison with the data. The figure shows how the two components of the radio LP (HEGs and LEGs) can be explained 
by a population of high-accretion rate SMBHs and one of low-accretion rate SMBHs, different distributions of spins. 



SO that any contribution from lower mass SMBHs must be 
small. 

The QSO models do not fit the entire HEG radio LF, 
but an important positive result is that they do not overpre- 
dict the density of lower-radio luminosity HEGs. Indeed, an 
extra population is probably needed to explain the higher 
space density of HEGs, and these are most likely SMBHs 



with lower masses and moderate or low accretion rates, 
for example Seyfert galaxies and lo w-luminosity AGN (e.g. 
iNagar et al.ll2005l : iFilho et al.|[2006h . 

We can also estimate the distribution for all SMBHs, 
by combining the spin distributions of ADAFs and QSOs 
weighed by their space densities: 
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Figure 5. The entire observed z=0 radio luminosity function down to 10^^ W Hz~^ sr~^ compared to our models with their best-fitting 
parameters. This figure is similar to Figure]?] but it shows a wider range in radio luminosity density. The parameters have only been fit in 
the region 10^^ W Hz~^ sr~^, so we are now comparing the predictions from our best-fitting models to the observed radio LF. We can 
see that the ADAF models provide a very good match to the LEGs all the way down to luminosity densities of 10'^^ W Hz~^ sr~^. There 
is no need for a further component to explain the LEGs, so our model predicts that the LEGs down to 10^^ W Hz~^ sr~^ are all powered 
by SMBHs with m,^ 10* Mq, with a negligible contribution from lower mass SMBHs. The models do not provide such good fits for the 
HEGs, indeed an extra component is required to explain the high space density of HEGs at luminosity densities < 10^'^ W Hz~^ sr~^. 
The most likely explanation is that SMBHs with masses below 10* Mq and moderate accret ion rates dominat e this region of th e HEG 
radio LF (these sources are typically Seyfert galaxies and low-luminosity AGN, LLAGN, e.g. iNagar et aLll2005l: iFilho et al.ll2006l ) . 



PsMBH(a) = 
where: 



t"QSO-PQSo(ffl) + t»ADAF-PADAF(a) 
WQSO + lUADAF 



WQSO = (t>QSo{Lv) 



(31) 



(32) 



WADAF = 4>AI>Af{Lv) 



(33) 
(34) 



The expectation value for this distribution is given by: 
{a)sMBHs = / aPsMBH(a)da, (35) 
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Figure 6. Spin distributions at z=0 for all SMBHs resulting from fitting the two Gaussian model to the local radio luminosity function. 
For each of the jet-efficiencies, the distribution shown is the combination of the ADAF-LEG and the QSO-HEG components quoted in 
Table[3] weighted by their relative space densities. Since the observed space density of LEGs is higher than that of HEGs, the distributions 
shown here are dominated by the ADAF-LEG component, with a small contribution from the QSO-HEG component. Note that the M05 
distribution has a component peaking at a~0.16, with a very small amplitude. 



The inferred distribution of spins for all SMBHs are 
shown in Figure[6l where we have combined the distributions 
for QSOs and ADAFs, and weighted them by the relative 
modelled space densities of QSOs and ADAFs respectively. 
All the distributions are bimodal, although the relative im- 
portance of the two components varies greatly between dif- 
ferent jet efficiencies (e.g. contrast HK06 and M05). Note 
that we weigh by the model distributions (jiQSoiLv) and 
4iadaf{L^), so that the expectation value will be most reli- 



able for jet efficiencies that can reproduce the observations 
better. 

For any given jet efficiency, the range of expectation 
values for the mean spin of all SMBHs is similar for models A 
and B as it is for model C (see Table |4] and Figure [6l). 

In Section [S] we discuss the implications of such a bi- 
modal distribution, and compare it to theoretical predictions 
from accretion and mergers of SMBHs. 

Interestingly, regarding the radio-loudness of quasars, 
our results suggest that there should be something resem- 
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bling a dichotomy in radio- loudness (see Figure|6]). There has 
been some debate r egarding the existence of a dichoto my in 
radio-loudness (e.g. Ilvezicl l2002l : ICirasuolo et al.]l2003h . but 
we note here that the dichotomy in spins will translate into 
a dichotomy in jet efficiency, which does not necessarily map 
linearly into a dichotomy in the ratio of radio-to-optical lu- 
minosity density (often used to define radio-loudness). The 
uncertainties in converting from jet power to radio luminos- 
ity density, the scatter in bolometric corrections and selec- 
tion effects can all blur out such a dichotomy. 



5.2 A prediction for the z=l radio LF 

If we assume that the individual spin distributions for 
ADAFs and QSOs do not change significantly up to red- 
shift of 2=1, then we can compare our predictions to the 
observed radio LF at this redshift. 

We can use the observed X-ray LF, which tells us how 
the space density of SMBHs with high accretion rates in- 
creases with z. However, we need to make an extra assump- 
tion about the ADAFs, represented by the BHMF. If we 
assume that the space density of SMBHs accreting at very 
low rates has not changed significantly between z=l and 0, 
then we can use the local BHMF, with the same distribu- 
tion of Eddington ratios (this is supported observationally 
by the relatively const ant evolution of AG N with moderate 
radio luminosities, see ISmolcic et al]|2009l '). 

The larger volumes probed at z=l, together with 
the strong evolutio n of the radio LF of FRII sources 
(|Willott et al.l 1200 ll ). mean that we are able to ob- 
serve radio sources up to the highest radio luminosities, 
~102* W Hz-i sr"^ at 1.4 GHz. Such sources have are typ- 
ically FRIIs and are rarely seen in the local radio LF due to 
the smaller local volume. 

Figure [7| shows the resulting radio LFs from the QSOs 
and ADAFs, compared to the observed radio LF, for the 
same range of radio luminosities that we have observed. We 
can see, that unlike at z=Q, it is the QSOs that dominate 
the radio LF. The solid line in Figure [7] shows the observed 
radio LF at z=l, combin ing the high radio luminosity LF 
Willott et al.l (12001 ) and the low radio luminosity LF 



Smolcic et all (|200 



I 



from 
from 

Our model overpredicts the radio LF in some regions, 
particularly around ~10^^ W Hz~^ sr~^. We note how- 
ever, that this is an ill-constrai ned region of the rad io LF: too 
faint for the samples used bv lWillott et al.l (12001^. but too 
rare in space density for the area covered bv'Smol cic et al] 
(|2009l ). Indeed the apparent 'wiggle' is a conseque nce of the 
analyt ic function fitted to the radio LF by IWillott et al.l 
l|200ll ). 

The good agreement over so many decades of radio lu- 
minosity is very encouraging. Looking at the results for the 
different jet efficiencies, we can make the following predic- 
tions for the 2=1 radio LF: 

At radio luminosity densities J^IO^® W Hz~^ sr~^, the 
composition of the z=l radio LF should comprise ~100% 
HEGs. This prediction is in very good agreement with ob- 
servations (Fernandes et al. 2011, in prep.). In the range 



10^ 



10 W Hz sr , at the higher luminosities we ex- 



pect the HEGs to dominate, with LEGs making up only 
~10% of the total population. We then expect the fraction 
of LEGs to increase rapidly with decreasing luminosity and 



eventually dominate. The dominance of the HEG popula- 
tion, unlike in the local radio LF, is due to the strong in- 
crease between z=Q and 1 of the comoving accretion rate 
onto SMBHs, whi ch is reflected in the s trong evolution in 
the X-ray LF (e.g. [Silverman et al.lliooil '). 



5.3 The comoving accretion rate of radio sources 

Applying Equations 1311 1341 and [35] to the accretion rate, 
rather than the spin, we can estimate the expectation value 
for the accretion rate onto ADAFs and onto QSOs. 

Figure [9] shows the evolution of the comoving accretion 
rate (left panel), and the fraction of the total accretion rate 
onto ADAFs and onto QSOs (right panel) as a function of 
redshift. At z=Q, approximately half of the accretion onto 
the most massive black holes is occurring via high accretion 
rates (in the form of QSOs), and the other half via low ac- 
cretion rates (ADAFs). The fraction of accretion ocurring in 
high-accretion rate objects increases rapidly with redshift. 
The right panel of Figure [51 shows our estimate, although 
again if the mass function of SMBH decreases significantly, 
the curve for ADAF will overpredict the accretion rate onto 
low-accretion rate sources. 

Our left pa nel of Figure [9] is very similar to Figure 4 of 
[Merloni fc Hein z (2008). For both works, the accretion rate 
from QSOs is derived from the X-ray LF, so unsurprisingly 
the curves agree very well. 

The difference lies in the density of accretion onto 
ADAFs, where we have assumed a distribution of Eddington 
ratios for the local SMBH mass function, and assumed this 
to stay constant with redshift. We have used de inferred dis- 
tribution of spins to pro duce a distribution of jet efficiencies. 
[Merloni fc Heina (|2008[ '). on the other hand, assume a fixed 
fiducial jet efficiency and derive the accretion rate from the 
radio LF. 

Similarly, [Shankar et al.[ (|2008[ ) use a sample of 
optically- and radio-selected AGN to derive the mean jet 
efficiency, together with the scatter. The resulting jet effi- 
ciency for the sample is ~0.01. The sample is actually the 
same one used in our Section [3] so that the same caveats 
mentioned in Section [3.21 should be bore in mind, since the 
sample includes only objects which are bright enough to be 
detected both in the optical and radio surveys. 



6 THE EVOLUTION OF THE COSMIC SPIN 

In this Section we use our results from Sections [J] and [S] to 
infer our best estimate for the cosmic spin history of the 
most massive black holes. Comparison of Figures [6] and [8] 
shows that as the redshift increases, the fraction of SMBHs 
with high spin decreases. 

We stress that Figure [8] assumes negligible evolution in 
the black hole mass function. This is probably an acceptable 
approximation up to z^l, but at higher redshifts we expect 
the space density of the most massive black holes to have 
decreased significantly. At higher redshifts, the distribution 
of SMBHs will closely match that of the HEGs alone, and 
most of the probability density will lie at low spins. 

Before we continue, we warn of two important points. 
First, it is not the intention of this paper to judge the dif- 
ferent models of jet formation around SMBHs. Second, as 
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Figure 7. The observed radio LF at 2=1 (solid line, IWillott et alll200ll : ISmolcic et al.ll200g^ compared to the predicted radio LF from 
our model. We are assuming that the spin distributions for high- and low-accretion rate SMBHs (QSOs and ADAFs) at z=l are the 
same as those at 2=0. The space density of QSOs is derived using the X-ray LF, at 2=1, while the space density of ADAFs is derived 
assuming the present day black hole mass function to be already in place at 2=1. The spin distributions of the QSOs for the BB09 and 
HK06 models lack a high-spin component (see Figure IBJ, hence the sharp drop at high radio luminosities. The overall agreement is very 
good, except around ~10^'' W Hz~^ sr~^, although we note that this r egion is not we ll constrained observationally, and the 'wiggle' 
is a consequence of the functional form used to describe the radio LF (Willo tt et al.|[200j) . Our model predicts that the radio LF at 2=1 
is dominated by high-accretion rate SMBHs, which should be observed as sources with HEG spectra. 



discussed earlier, the quality of the fit is partly determined 
by our choice of / factor when converting jet power to ra- 
dio luminosity, and a lower value would cause the models 
that currently look best to overpredict the radio LF, so that 
some of the less-favoured models would then provide better 
fits (see discussion in Section [5. ip . 

Hence we choose to estimate the mean cosmic spin his- 



tory using two sets of efficiencies. The first uses all the mod- 
els independently of whether they can reproduce the radio- 
loudness of quasars or the local radio LF. The second will 
use only the four efficiencies that show some success at re- 
producing the observations, namely TlOA, BB09, N07 and 
HK06. We label these as the 'best fitting set'. 

The cosmic spin of SMBHs evolves from distributions 
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Figure 8. Spin distributions at z=l for all SMBHs. These have been obtained in an analogous way to Figure|6l only using the X-ray 
LF at 2=1. The vertical dotted lines mark once again the mean value. We have assumed the ADAFs not to evolve significantly up to 
z=l, which means we are likely to overestimate their contribution. Comparing this figure with Figure [B] we can see that the fraction of 
SMBHs with a high spin decreases with increasing redshift. 



dominated by low spins at high redshifts, to bimodal distri- 
butions with comparable high-spin and low-spin components 
at low redshifts. It is the fraction of highly spinning black 
holes that increases with cosmic time. To quantify this, we 
consider the fraction of SMBHs with a>0.5. The fractions, 
or probabilities of a^O.5, are shown for each jet efficiency in 
Table [5l together with the mean fractions of all the sources, 
and of the best fitting set only. If we consider all the efficien- 
cies together, the fraction of SMBHs with increases 
from 0.34±0.29 at z^l to 0.45±0.33 at z =0. However, we 
believe a more accurate result is given by the best fitting 



set of efficiencies, TlOA, BB09, N07 and HK06: in this case 
the fraction of SMBHs with is 0.16±0.03 at z=l and 

0.24±0.09 at z=0. The uncertainties are due to the varia- 
tions of the fraction between different efficiencies (see Ta- 
ble [SJ. At redshifts >1 we expect the fraction of SMBHs 
with to be even lower than the quoted values. 

Given the shape of the distribution functions, the most 
informative way of describing the spin history is to show 
the distributions and determine the fraction with high spin. 
A mean value does not reflect the bimodality, however, for 
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Figure 9. Expectation value for tlie comoving accretion rate from QSOs and ADAFs (left panel) and relative fraction of the accretion 
rate (right panel), derived in an analogous way to Equations 1311 l34l and Equation |35] but applied to rh rather than d. 
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P(a^ 0.5) 
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0.14 


0.13 
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0.94 


0.78 


BB09 


0.33 


0.19 


N07 


0.20 


0.13 


HK06 


0.30 


0.17 


M05 


0.80 


0.62 


all 


0.45±0.33 


0.34±0.29 


best fitting 


0.24±0.09 


0.16±0.03 



Table 5. Fraction of SMBHs with a>0.5 at z=0 and 1, for each 
of the six jet efficiencies. At the bottom we show the mean values 
for all the efficiencies and for the best fitting set only. 

completeness we have calculated the expectation value and 
its evolution, and we describe this in Sections 16. II and 16.21 

6.1 The expectation value of spin with redshift 

Figure [To] shows the evolution of {a{z)) (solid line) together 
with the standard deviation (dashed lines). The standard 
deviation is measured from the square root of the variance. 
The large value of the variance is caused by the bimodal 
nature of the spin distribution at each redshift (as shown by 
Figures |S] and [S| . 

Due to the fact that the spin distributions are bimodal, 
the expectation value should be taken simply as a fiducial 
value, and the exact value is clearly uncertain. However, this 
expectation value evolves with redshift, and we see that it 
was lower at z=l and has slightly increased to a higher value 
in the present-day Universe (see Figure [TD}. In Section [5] we 
discuss the physical implications of this evolution. 

6.2 A best estimate of the evolution of the mean 
spin 

Figure [To] shows a range of different spin histories, due to the 
differences in the jet efficiencies from Figure [2] We wish to 



use all the available information, in the form of the different 
jet efficiencies, for a best estimate of the spin history. We 
consider both the entire set of jet efficiencies, as well as 
considering the best-fitting set alone. 

To obtain the best estimate of the spin history, we 
take the mean of all the results as a function of red- 
shift. However, the individual results are not always com- 
patible. The A'^ individual results di{z) are compatible if 

[{ai{z) — {a{z))) joif is ^(N — 1), wher e (a(z')) is the 
mean from all the individual results (see e.g. |Press|[l997l ). 

In order to combine incompatible results, we consider 
the probability that each measurement of ai{z') is correct, as 
well as the probability that it is not correct. A correct mea- 
surement is one where the quoted uncertainty is accurate, 
whereas an incorrect measurement is quoting uncertainties 
that are too small or might suffer from systematics. Incor- 
rect measurements have error bars that are too small, hence 
they are incompatible. 

We follow the method describe by jPressI (llQQTj ). and 
in particular their Equation 16, which gives the posterior 
probability for the value of (0(2)) given data: 

P((a(z))|data) oc ^ P(p) J|[pP« + (1 - p)Pn]dp (36) 

i 

where Pci and Pu are the probability distribution functions 
for correct and incorrect estimates, and the product J^, is 
over the rii different estimates. 

The first term in the square bracket is the product of 
the probability p that measurement i is correct times the 
probability distribution function for that measurement (us- 
ing the standard deviation of the correct measurement, aci)- 
The second term is the product of the probability 1-p that 
the measurement is incorrect. The probability distribution 
function Pu therefore requires a standard deviation which 
is large enough to make the measurement correct. We label 
this as a\i. The integral over dp means we marginalise over 
both cases. 

We approximate the probability distribution functions 
for (0(2)) as gaussians. In the case of Pci , we assume the 
standard deviation aci to given by the errors shown on Fig- 
ure [10] 
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Figure 10. Expectation value for the spin of SMBHs as a function of rcdsiiift (solid line), as determined by combining the spin 
distributions for ADAFs and QSOs, using the X-ray LF and black hole mass function to weigh the relative contributions. The dashed 
dotted lines show the values of ±1 standard deviations, where the standard deviation is the square root of the variance determined from 
the weighted spin distribution. In the TlOB, BB09, N07 and HK06, panels the -Icr contour fall bellow a=0. In the case of M05, the +la 
contour is higher than a=l. 



We now consider two extreme cases: 

1. There is the possibility that some of the estimates 
of {a{z)) are indeed incorrect. If, however, incorrect esti- 
mates were assigned a larger uncertainty, they would be- 
come correct. How large this new uncertainty should be is 
somewhat arbitrary, and might vary from model to model. 
Hence, we assign the same very large new uncertainty to all 
estimates, aii=ai —1, and within this uncertainty all our es- 
timates agree. We also assume a flat prior for the probability 



that measurements are correct, in the range O^p^l, so that 
P(p) = l. Equation [36] then simplifies to: 

P((a(z))|data) = 

^ + (37) 

2(2^)-n,<^c> 2(27r)- 

2. The measurements are all correct, since their uncer- 
tainties are large, they are formally consistent with each 
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Figure 11. Best estimates of tiie evolution of the mean cosmic spin of SMBHs with redshift, derived using Equation 1361 The left panel 
combines all the results of Figure [TOl and since these are not compatible, it considers the probability that each independent estimate is 
correct as well as the probability that it is incorrect, and marginalises over these (see Equation I37I I. The definition of compatible results 
is given in Section l6.2l The right panel uses the results from TlOA, BB09, N07 and HK06 only (the 'best fitting set'), and since these are 
compatible, it considers all the independent estimates to be correct, and is simply the product of their individual posterior probability 
distribution functions (Equation 138 II . In this case a more modest spin evolution is seen between z=l and 0., and note that the -la and 
-2(T contours fall mostly below the a=0. Both estimates assume that the space density of low-accretion rate SMBHs, ADAFS described 
by the local black hole mass function, is constant with z. Since the ADAFs tend to have higher spin values, this approximation is likely 
to overestimate the mean spin at the high-redshift end (z<;l). 



other. We therefore assign a delta-function prior at 'cer- 
tainty', P{p)—5{p—l) and Equation 1361 simplifies instead to 
the standard result: 

P((a(2))ldata) = (38) 

(2^)^11; ^c. 

Figure [TT] shows the resulting evolution of {a{z)} under 
these two assumptions. When combining all the results of 
Figure [TOl some of the results are incompatible, so that we 
use case 1 (Equation I37|) . The typical spin at z—1 is a~0.35 
and increases to ~0.75 at z—0. The right panel of Figure [TT] 
shows the resulting spin history when combining the TlOA, 
BB09, N07 and HK06 results only, using case 2 and Equa- 
tion 1381 In this case a more modest increase is seen, from 
a~0.25 at z—1 to ~0.35 at z—Q. The right panel is derived 
from our best fitting set, and represents our best estimate 
of the mean spin history. 

We stress that at all redshifts, we are assuming the pop- 
ulation of ADAFs is well represented by the z=0 SMBH 
mass function convolved with the distribution of Eddington 
ratios discussed in Section [4.21 This is probably a good ap- 
proximation up to z~l, since most of the mass of SMBHs 
with m.^lO** M0 was accreted at higher redshifts. The 
SMBIf mass function is probably similar at z~l and z~0, 
but at z^3 it will be significantly different. The space den- 
sity of high-mass SMBHs sould be much lower. Hence our 
inferred spin history contains too high a contribution from 
ADAFs at high redshifts, so that we probably overestimate 
the spin at z^l. 

Finally, we note that although the evolution of the ex- 
pectation value, {a{z)), appears weak, this is not necessarily 
the best way to illustrate the evolution in spin for such bi- 
modal distributions. The evolution in spin is stronger than 
what {a{z)} would suggest: at high redshift most of the 
SMBHs have the spin distribution describing the QSOs so 



that most sources have very low spins, with a negligible frac- 
tion having high spins. At low redshift, most of the SMBHs 
are described by the ADAF spin distribution, meaning that 
around half of the sources have high spins. 

It is important to note that the evolution in the spin is 
driven by the disappearance of high-accretion rate, low-spin 
SMBHs, whose space density decreases enough that the ap- 
proximately constant population of low-accretion rate, high- 
spin objects becomes dominant. 



7 COMPARISON TO OTHER INFERENCES 
ON BLACK HOLE SPIN 

The spin of black holes in AGN have been estimated in dif- 
ferent ways, and our results do not necessarily agree with 
those of other gro ups. For example, lElvis et al.l (|2002l ) and 
IWang et al.l (|2006l ) have claimed that most SMBHs must be 
spinning rapidly, while more recently iFender et al.l (|2010l ) 
have found no correlation between black hole spin and jet 
power for X-ray binaries. We therefore compare our results 
to previous work and discuss the possible reasons for dis- 
crepancies. 



7.1 The mean radiative efflciency 

We begin by comparing to estimates of the spin from the 
mean radiative efficiency. Comparing the total energy ra- 
diated by AGN over cosmic time with the mass density of 
SMBHs in the local Universe can be used to estimate the 
mean radiative efficiency (Soltan 1982). The mean radiative 
efficiency can be converted into a black hole spin using the 
NT73 model, but as is shown in Figure [T] the radiative effi- 
ciency from simulations can differ from the NT73 by factors 
of ~2, resulting in large uncertainties. 
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lElvis et alj (|2002l ) used the intensity of the X-ray back- 
ground and assumed the bulk of the AGN causing it to be 
at z=2. The authors inferred {e)^0.15, suggesting the mean 
spin of SMBIfs is a~0.9. The bulk of the sources producing 
the X-ray backround, however, are at lower redshifts than 
their assumed value of z=2. A more appropriat e characteris- 
tic redshift would have been z=0. 7-1.0 (see e.g. lBarger et alJ 
l2003l : [Ueda et al] |200vj). which would bring their mean effi- 
ciency down to values of {e)<;0. 08-0.1. 

Ano ther high va l ue of the radiative efficiency was in- 
ferred bv lWang et all (|2006l ) using a different method. These 
authors compared the luminosity function of quasars to the 
mass function of quasars in redshift bins, and estimate (e) to 
be in the range 0.20-0.35, suggesting the SMBHs are spin- 
ning rapidly. 

In order to derive a radiative efficiency, IWang et al.l 
()2006l ) equate the integral of the luminosity function 
(/ iboi<^(iboi)dI/boi) to the integral of the mass-density 
function of quasars {J m,(f)(m,)dm,). However, these au- 
thors derive the mass functions from the broad-line widths 
of the quasars in a magnitude-limited sample. Only quasars 
brighter than the limiting magnitude (_B-band magnitude of 
19.4) can have their masses measured and included. Since 
quasars have a range of Eddington ratios (e.g. McLure & 
Dunlop 2004), the correspondence between bolometric lu- 
minosity and black hole m ass is not one-to-on e. While the 
luminosity integral used by IWang et al.l (|2006h is complete 
above their flux limit, the mass integral is not. This is be- 
cause not all quasars with mass above their nominal mass 
lower limit are optically bright enough to be detected in 
their survey so that their mass function is incomplete. The 
fact that the luminosity integral is complete but the mass 
integral is incomplete means that they will overestimate the 
amount of energy that was radiated away. We therefore con- 
sider their estimate of the radiative efficiency, and hence 
spin, as an upper limit. 

Most of the recent estimates of the mean efficiency 
from the Soltan argument typically find values of (e) in 
the range 0.06-0.1, with a typical v a lue of (£)~0.07 (e.g . 
Marconi et all l2004l: ICao fc Lil l2008l : IShankar et aD l2009l : 



Martmez-Sansigre fc TavloJ 2009[ see their Section 5 for 



a discussion of recent results). 

We can compare all the previously-mentioned estimates 
to our results by computing the luminosity-weighted mean 
spin: 



(a)i 



/ (a(z))best{iboi)d2; 



/ {Lboi)dz 

/ [/ a(2)P((a(z))|data)da] [/ CxL^^Lx, z)dLx] dz 
J J CxLx(t){Lx, z)dLxdz 



(39) 



For the weighting function in Eauation l39l we use the X- 
ray LF in the range OsJz^S, and due to the strong evolution 
of the X-ray LF, practically all of the weight lies with the 
high redshift values of {a{z)}, since the bulk of the luminosity 
was emitted at z^l. 

P({a(2))|data) is given by the left panel of Figure [TTl 
which then yields a value of {a)L=0.40. The resultant mean 
radiative efficiency (e)L is therefore e(0.40)=0.075. However, 
when using the spin history from the right panel in Figure fTTI 
(the best fitting set), we infer {a)L=0.28, or e=0.068. The 



mean radiative efficiency derived from our spin history is 
therefore in excellent agreement with the estimates from the 
literature. 



7.2 Lack of correlation betvireen estimated spins 
and jet powers in galactic black holes 

Comparing the spin measurements of GBHs (ga l actic black 
holes) with their radio luminosities. iFender et all (|20ld ) con- 
clude that there is evidence for an additional parameter 
which controls the amount of power in jets, yet this param- 
eter does not correlate with the spins estimated from fitting 
the accretion disc or the reflection edge. 

In many a spects, AGN behave a s scaled up versions 
of GBHs (e.g. liviarscher et al.l l2002l: iMerloni et all l2003l : 



iFalcke et aLll2004l : iMcHardv et al.ll2006l '). The lack of corre- 
lation between 'radio loudness' and spin in GBHs is therefore 
a concern for our work. 

It is encouraging that when plotting a measure of bolo- 
metric luminosity versus a measure of jet power, both GBHs 
and SMBHs show correlations spanning several decades, but 
also a scatter of 2 or 3 decades in r adio luminosity (e.g. 
ISikora et al] l2007l : iFender et al]|2010l ). In this paper, the 
scatter has been interpreted as arising mainly due to varia- 
tions in black hole s pin (see also lSikora et al.ll2007l ). however 
IFender et all (|2010l ) find no correlation with published spin 
estimates for the GBHs. IFender et all l|2010D conclude that 
the lack of correlation could be due to spin having no effect 
on the observed jet power, or alternatively due to either the 
spin measurements or the measurements of jet power being 
inaccurate. 

It is possible that at low spins, the power from the disc 
jet is larger than the power from the black hole jet (e.g. the 
HK06, N07 and BB09 jet efficiencies in Figured]), compli- 
cating the interpretation. However, if spin has an observable 
effect on the jet power, then one would expect some correla- 
tion between spin and jet power. If the disc jet were always 
more important than the jet driven by the spin of the black 
hole, then the spin would be irrelevant when explaining ob- 
servations of jets. 

Subsection 17.2. 1 1 discusses the uncertainties in estimat- 
ing the jet power from core radio luminosities. Subsec- 
tion 17.2.21 discusses the uncertainties in the estimates of 
black hole spin from fitting the accretion disc luminosity 
or the iron refiection lines, and Section 17.2.31 discusses the 
observed lack of correlation. 



7.2.1 Estimated jet powers 

IFender et all l|2010l ) measure the jet powers from the core 
luminosities, as is usual in the GBH community. How- 
ever, for AGN the core luminosity is totally incapable 
of predicting the total radio luminosity (see Figure 7 
of iKlockner et al.l |2009| ). and it is the total radio lu- 
minosity that is required to trace the jet power (e.g. 
iRawlings & Saunders 199l'; ' Willott et al]| 19991 : iBirzan et all 
2004; Cavagnolo et al... 20101 

We note the cases of Cygnus X-1 and S26 (a micro- 
quasar with extended radio structure reminiscent of an FRII 
radio galaxy) where the extended lobe radio lum inosity is ~2 
decades brighter than expected from the cores (|Gallo et al.l 
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l2005l : ISoriaet al.ll2010l ). which in the case of S26 is not even 
detected. It seems that in GBHs, the core luminosity might 
not be a good tracer o f the total jet power either (see also 
e.g. iKaiser et al.ll2004l . for the case of GRS 1915+105). 

To complicate things furth er, the discovery of a one- 
sided bow-shock in Cygnus X-1 (|Gallo et al.ll2005l ). suggests 
the environments of the GBHs can be much less homoge- 
neous than those of AGN, on the respective scales of their 
jets. 

The extended jets will not be detectable unless they 
are seen to do work on the interstellar medium (as is only 
the case in the Northern lobe of Cygnus X-1). This suggests 
accurate determination of the total jet power in GBHs might 
suffer from even larger uncertainties than in the AGN case 
(Section EH). 



7.2.2 Estimated spins 

The spin estimates used bv lFender et al.l ()2010l ) are obtained 
either from modelling the thermal continuum from the ac- 
cretion disc ( Zhang et al.,1997,) or from the iron fluorescence 
reflection lines ( Fabian et al.l 19891 ) . 

The spin measurements from the iron reflection (K-a 
and L) lines assume that the red wing of the line is due 
to general-relativistic effects in the vicinity of a black hole. 
From the line width, the characteristic radius where most of 
the reflection is occurring is estimated. 

The disc-fitting methods model the spectrum of the 
source as continuum thermal emission originating in an ac- 
cretion disc. The free parameters include the mass of the 
black hole and the spin, which determines the truncation 
radius, the accretion rate as well as the tilt angle of the 
accretion disc. 

The two methods rely on identifying a characteristic 
radius. This radius can be where the continuum emission 
truncates or where the bulk of the reflection originates. If 
the observed characteristic radius corresponds closely to the 
ISCO, which depends on the spin parameter in the Kerr 
metric, then the characteristic radius yields an estimate of 
the spin. 

However, if the characteristic radius does not corre- 
spond to the ISCO, then the resulting spin estimate will be 
incorrect. One concern is therefore whether significant radi- 
ation or reflection can occur at radii smaller than the ISCO, 
so that assigning the characteristic radius to the ISCO un- 
derestimates the ISCO and therefore overestimates the spin. 

Several authors have studied this possibility, with vary- 
ing conclusions, a n d we summarise here recent progress. 
[Ab ramowicz et all (|20ld ) discuss how different methods 
(disc-fitting, iron lines) will probe different characteristic 
radii. They argue that for low Eddington ratios, A<0.3, im- 
plying thin discs, the approximation that these radiation or 
refiection edges will correspond to the ISCO is valid, but that 
this approximation will break down for the hig her Eddington 
ratios that produce thicker discs. The results of lSteiner et al.l 
l|2010ll for LMC X-3 support this: analysing data that span 
26 years, and using only the data for periods during which 
the Eddington ratio is 0.05s;A^0.3, they find an inner disk 

radius that is constant to within 4- 6%. 

The simulations b y IShafee et al.l l|2008D . 



iRevnolds fc FabianI (120081 ) and lPenna et all l|2010l ) also sup 



port the view that for thin discs (/i/r;$0.1) the deviations 

from NT are small. 

The simulations of iBeckwith et^ l|2008bl ) and 
iNoble et al.l ([2009 ) suggest that deviations from the 
NT73 model are pa rticularly important at low-spins. 
IBeckwith et al.l (|2008bl ) also discuss the possible dependence 
of this on h/r, and note that height in their simulations is 
~0.06-0.2. 

For modera tely thick discs (/i/r~0.2), the simulations 
bv lFragild (|2009| ) found that when the tilt angle was zero, the 
dependence of the characteristic radius with spin resembled 
closely that of the NT73 model. However, with tilt angles 
of only 15 degrees, the characteristic radius is constant with 
spin, and any attempt to infer the spin from the charac- 
terist ic radius will severely underestimate the spin (|Fragilel 
I2OO9I ). 

Another potential problem can occur if the accretion 
disc truncates at a larger radius than the ISCO, for example 
if the accretion rate is too low to sustain an optically-thin 
disc in the innermost region, and the accretion flow becomes 
an ADAF instead. In this case, the radius of the ISCO will be 
overestimated, meaning the inferred value of the spin will be 
a lower limit. Indeed, in Cygnus X-1, the fits to the thermal 
continuum and the K-a line suggest the tru ncation radius is 
at a distance ^2-3 x larger t han the ISCO (|Gierliriski et al.l 
I1999I : iMakishima et al.l l2008'fl 

For one of the s ource s with a measured high spin 
(a=0.94 iMiller et al.l |2009| ). a subsequent analysis by 
iDone fc Diaz Trigol (|2010l )~ has found that the extremely 
broad line is artefact due to the brightness of the source 
causing 'Pile-up' of photon^. The authors find that after 
taking this effect into account the line is narrower than pre- 
viously determined, corresponding to a characteristic radius 
that does not require high spin and might even imply a trun- 
cated disc. 

Indeed, we note that spin measurements of the same 
GBH from fitting the disc continuum and from the fitting 
the reflection components do not alw ays agree within th e 
quoted uncertainties (see Table 1 of iFender et al.|[201ol ). 
Some GBHs have had their spin estimated both from from 
disc-fitting and from iron lines, yet the two estimates do not 
always agree within their formal uncertainties. 

From disc-fitting alone, the reported measurements for 
GRS1915-I-105 from dis c-fitting range from t o 0.998 (see 
refere nces in Tender et al.ii2010l ). More recentlv. lBlum et al.l 
(|2009l) have estimated the spin again from disc-fitting and 
from the K-a reflectio n line, and obtain s pin values of a=0.56 
and 0.98 respectively. iBlum et al] (|2009l ) use observations of 
GRS1915-(-105 during a low-hard state only, with A<0.3, 
and argue that some of the discrepancy with other groups 
could be partly due to observing at different states (e.g. 
iMiddleton et al.ll2006l . observed the source when A>0.3). 

Finally, an additional complication arises from the pres- 



^ For another example where the t runcation radius mi g ht lar ger, 
see the study of XTE J1650-500 bv lDone fc Gierlinskil 1I2OO6I) . 
^ Pile-up occurs when more than one photon reaches the CCD 
in between two read outs. The local excess in charge will appear 
larger, hence a higher energy will be erroneously assigned to the 
photon. However, if the excess is above a cut-off t hresho l d, the 
fl ux from photons will not be counted at all. See e.g. lBaUe3 lll999l ) 
or lDone fc Diaz Tried 1I2OIC1I) for more details. 
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ence of absorbing material along the line of sight and 
near the black hole. Photoelectric absorption from gas with 
a range of ionisation parameters can lead to a spectral 
shape similar to a red wing, and some authors suggests 
that the broad reflection lines can be fitted without re- 



Miller et al. 


2008; 


Turner & MilleT 20091: Mivakawa et all 


20091: Yaaoob et alj 


2010l~). This is an area of intense debate. 



where other authors suggest that rela tivistic blurring is in- 
deed required to explain the lines (e.g. [Reynolds et al.|[2009l : 
iFabian et ai]|2009l : IZoghbi et allboioh . 



7.2.3 



Lack of correlation between spins and jet powers 



To summarise, while measurements from the reflection lines 
of disc fitting offer an exciting and promising method to 
determine black hole spin, the accuracy of these measure- 
ments is still being debated. Our primary concern is that the 
characteristic radius inferred by these measurements might 
not always correspond to the ISCO. These concerns apply 
to fitting both the continuum-fitting and the refiection line 
methods. 

This is not a criticism on our part, since the spin es- 
timated from rj in our work is also hugely uncertain and 
Figure [2] shows the difi'erent jet efficiencies we have used: 
it is clear that they vary enormously. Converting between 
jet powers and radio luminosities is also highly uncertain in 
GBHs as well as in AGN (see Sections ISTTI and ?L2l} . An 
additional consideration to bear in mind is that measures of 
AGN jet power from the total radio luminosities are actually 
derived from the total energy deposited over typically ~f0^ 
years, so that they are only time-averaged values. 

Given the uncertainties we do not think that the ob- 
served lack of correlation be t ween spins and radio luminosi- 
ties found by I Fender et al.l (|2010l ) is at odds with our as- 
sumption that spin is a monotonically-increasing function 
of spin. The uncertainties in the spin estimates and the es- 
timates of the jet power are still too large. We do not think 
that the current evidence is strong enough to conclusively 
rule out the spin paradigm. 

We stress that there is strong evidence amongst GBHs 
for the je t power to correlate with the state of the accretion 
flow (e.g. lFender et all 12004 iDone et al.ll2007l '). We are not 
arguing here that such effects are not important in AGN. 

Amongst AGN, however, huge differences in (time- 
averaged) jet power are also observed with no noticeable dif- 
ference in the accretion flow (reflected, for exa mple, by the 
shap e of the SEP or the Eddington ration, e.g. lElvis et al.l 
1 1994 ISikora et all 120071 ). It is this difference that requires 
an additional hidden variable such as spin. 

Finally, we note also that winds with similar escape 
velocit ies have been detecte d in radio-loud and radio-quiet 
AGN (|Tombesi et al.ll2010al . 2010b). Their results are con- 
sistent with both populations having similar underlying dis- 
tribution, but the radio-loud sample is much smaller (5 ob- 
jects) than the radio-quiet sample (42), and unsurprisingly 
the latter shows a larger spread in velocities, with maximum 
values ~0.3c. These winds are likel y to be powered by ac - 
cretion close to the Eddington limit l|King fc Poundsll2003D , 
rather than jets. 



8 DISCUSSION 

We have used an assumption where the power available for 
the production of jets is a fraction of the accreted energy, 
rj. We have used the values for ri given by recent theoreti- 
cal and simulation work, summarised in Figure [21 This ap- 
proximation can explain the radio loudness of quasars as 
a dependence on the spin of the accreting SMBH as well 
as explaining several observed properties of radio-selected 
AGNs. 

In Section 2] we use three parametrisations for the 
spin distribution of high- and low-accretion rate SMBHs, 
and use the radio LFs of HEGs and LEGs together with 
a bayesian model-selection method to select the most ap- 
propriate parametrisation. A two-gaussian model is chosen, 
ahead of power-law distributions or single-gaussian models. 
This two-gaussian distribution of spins can explain reason- 
ably well the radio LF of HEGs and LEGs individually. 

It is interesting that both the high-accretion rate regime 
(in the form of the upper envelope for the quasars in Sec- 
tion [Sjl, as well as the low-accretion rate regime (the ADAF- 
LEGs in Section 3} require a mechanism to produce very 
high jet efficiencies, r/max^l. This seems to be independent 
of the accretion mode, and black hole spin provides one pos- 
sible mechanism for such high jet efficiencies, as shown by 
some of the models in Figure (2] 

For the two-gaussian distribution, the best-fit param- 
eters for the LEGs-ADAFs choses two gaussians of simi- 
lar amplitude, one with minimum spin and the other with 
maximum spin. This suggests that there are similar num- 
bers of high-spin and low-spin SMBHs amongst the SMBHs 
with low-accretion rates. The best-fit distribution for the 
HEGs-QSOs is also found to have a bimodal distribution, 
but this time with a very low amplitude for the high spin 
peak. Hence, most of the SMBHs with high-accretion rates 
have low spins. 

Since we infer a bimodal distribution for spins, which is 
heavily favoured over smooth power-law distributions, our 
model predicts a dichotomy in the radio- loudness of quasars. 
However, detecting such a dichotomy might prove difficult 
due to source-to-source variations in the conversion from jet 
power to radio luminosity density and from bolometric to 
monochromatic luminosity, as well as due to observational 
selection effects. 

We are able to make falsifiable predictions about the 
contribution of HEGs and LEGs to the radio LF. At z=0, 
we predict that the radio LF of LEGs is dominated by black 
holes with m.^10^ Mq down to radio luminosity densities 
of -10^^ W Hz"^ sr"^ at 1.4 GHz. At z=l and above a 
certain luminosity density, which we estimate to correspond 
to ~10^^ W Hz~^ sr"^ at 1.4 GHz, we expect aU the radio 
sources to be composed entirely of HEGs. Below this value, 
we expect a rapidly increasing number of LEGs. The ap- 
pearance of a large population of HEGs amongst the radio 
sources is due to the the strong increase of the comoving 
accretion rate onto SMBHs. 

As discussed in Section |6l the evolution in cosmic spin 
is most notable in the fraction of SMBHs with high spins. 
This fraction is low at high redshifts, and increases towards 
lower redshifts. The evolution in spin is really driven by 
the gradual disappearance of a high-accretion rate but low- 
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Figure 12. Distribution o f final black ho l e spin s from a Montecarlo simulation of 10^ SMBHs undergoing one major merger each. We 
use the analytic formula of iRezzoUa et al.l l|2008h . which requires the initial masses and spins of both black holes, as well as three angles 
describing the relative orientations with each other of the angular momentum vectors of each black hole and the total angular momentum 
vector (for a total of 7 free parameters). The three angles have fiat probability distribution functions between and 27r. The two initial 
spins are allowed any values between and 0.1, while the three angles any value between and 2it. The smaller mass is always fixed at 
a value of 1. In the left-hand panel, the larger mass is drawn from a flat distribution within the range 1-4. We are therefore representing 
what is normally considered as a major merger. In the right-hand panel, the larger mass is fixed at 1, representing equal mass mergers. 
The vertical dashed line represents the mean. The panels illustrate that high mean spin values, (a)=0.6-0.7 can be reached by single-event 
major mergers. Allowing the initial spins to have values O^a^Jl instead, with a flat distribution in this range, leads to much broader 
distributions of the final spins, but with essentially the same mean values. 



spin population, so that the low-accretion rate and high-spin 
population becomes gradually dominant. 

Combining the results from the different jet efficiencies, 
with the best-fitting parameters, and the evolution of the 
X-ray LF, we are able to estimate the expectation value 
for the spin as a function of redshift for black holes with 
10* M0: We find marginal evidence for evolution in 
the spin distributions. If we use all the efficiencies, we find 
the fraction of SMBHs with increases from 0.34±0.29 

at z—1 to 0.45±0.33 at z =0. If only the best fitting set 
is used, this fraction is found to be 0.16±0.03 at z—1 and 
0.24±0.09 at 2=0. 

The mean cosmic spin increases from a~0.35 at z^l to 
a~0.75 at z~0. However, when using the three best fitting 
jet efficiencies only, the evolution is found to be more mod- 
est: from a~0.25 to o~0.35 in the same redshift interval. 
This second set of values represents our best estimate. 

The luminosity-weighted spin comes out to be aL=0.40 
for all efficiencies and aL=0.28 for the best fitting set. This 
predicts the mean radiative efficiency to be (e) =0.075 (or 
(e)=0.068 when using the best fitting set). These quantities 
can be compared to independent constraints from the Soltan 
argument, and are in excellent agreement with them. 

Overall, we infer a very modest evolution of the mean 
spin, although given the bimodal distributions a mean value 
is not the best description of the evolution of the spin. In- 
deed, the spin distributions evolve from distributions dom- 
inated by a single low spin component at high redshift, to- 
wards bimodal distributions at low redshift. In addition, as 
stressed in Section [6l we are probably overestimating the 
high-spin contribution at Zftl- 

Perhaps counter-intuitively, the epoch of lowest black 
hole spin corresponds to the epoch of highest accretion. 
The evolution of {&{z)} inferred here, however, is consis- 
tent with the picture of 'chaotic accretion': matter from the 



host galaxy falls into the quasar from random directions and 
with different orientations of the angular momentum vector 
(see King et al.ll200^ '). 

Depending on the angle between the angular momen- 
tum vectors of the black hole and the accretion disc, these 
two will either align or counter- align. During approximately 
half of the time the accreted matter has the opp osite di- 
rection of angular momentum to the black hole (jKing et al.l 
I2OO5I : iLodato fc Pring"i3l2006l '). Co-rotating orbits are phys- 
ically closer to the event horizon of the hole than counter- 
rotating ones. Therefore co-rotating matter increases the 
angular momentum by a smaller amount than the decrease 
brought by counter-rotating matter. Chaotic acc retion is ex- 
pect to lead to SMBHs with low spins (a~0.1 iKing et al.l 
|200^ , and our result is in agreement with spun down AGNs 
at high redshift, when the mean accretion rate was higher. 

For accretion to significantly affect the spin o f a black 
hole, the increase in mass must be significant (|Bardeenl 
We approximate this as the condition 



which, if not met, implies no significant decrease in a due 
to chaotic accretion. 

The higher mean spin at lower redshifts suggests major 
mergers of SMBHs are driving the increase in spin. Galaxy 
major mergers are expected to lead to major mergers of the 
central SMBHs. During the coalescence of two black holes 
of similar mass, the orbital angular momentum will con- 
tribute significantly to the final angular momentum of the 
coalesced black hole. Hence, major galaxy mergers are ex- 
pected to lead t o central SMBHs with a moderately high 
spin value (e.g. iHughes fc BlandfordI l2003l : iRezzolla et alj 
I2OO8I I. Minor black hole mergers, on the other hand, tend 
to lead to low values of spin in the coalesced hole (see 
iHughes fc BlandfordI [2003 ). Values of (a)~0.7 can be rea- 
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Figure 13. Simulation of final black hole spins after 10® SMBHs 
undergoing a variable number of mergers. The number of merg- 
ers undergone is drawn from a Poisson funct ion with mean 0.7 
(following the results of iRobaina et al.l |2010|). The final spin i s 
calculated using the analytic formula of iRezzoUa et all l|2008l ) 
(see caption to Figure I12I I . We have assumed flat distributions 
for all the two black hole masses, the two black hole spins and 
the three angles. In this simulation, the black hole masses are 
assumed to be a constant fract ion 1.4 X 10~^ of the galaxy stel- 
lar masses l lHaring fc Rixlbooj) . These stellar masses are allowed 
to have masses in the range 5 X 10^" — 4 X 10^^ Mq, and are 
drawn at random from a Schechter function with a break at 
mass 1.4 X IQ-'^-'^ M(7), and low-mass end power law index of -1.2 
llCole et al.ll200ll) . The initial spins of the SMBHs are drawn from 
a flat probability distribution between a=0 — 0.1, and the spin of 
the coalesced SMBH is tracked (later incoming SMBHs are also 
assumed to have an initial spin in the range 0-0.1). The mergers 
are assumed to occur isotropically, and the flnal bl ack hole spins 
are ag ain inferred using the analytic formula of Rez zoUa et al.l 
1I2OO8I ). The distribution shows 3 obvious signals: the sources hav- 
ing undergone zero mergers are gathered around (a) =0-0.1, those 
having merged once peak around 0.65, and those having merged 
twice extend up to 0.85. The resulting expectation value of a=0.35 
is marked by the vertical dashed line. 



sonably achieved by SMBH major mergers, but values ~1 
are very hard to reach. 

In the low-redshift Universe, the most massive SMBHs, 
which dominate our estimate of {a{z)), are found in the 
most massive gala xies (with bulge masses ^ 10^^ M©, 
iHaring fc Rb3l2004l l. From mass-selected sam ples of galaxies 
and u sing the two-point correlation function. iRobaina et al.l 
l|2010l ) suggest that since z = 1.2, galaxies more massive 
than 1 X 10^^ M© have undergone, on average, 0.7 merg- 
ers involving two progenitor galaxies more massive than 
5 X 10^" Mq. 

To illustrate the effect of every SMBH undergoing one 
major merger we have carried out a simple simulation. Fig- 
ure [12] shows the distributions of final spins from two Mon- 
tecarlo realisations of 10® black hole single-event mergers. 
Once again, the amount of accretion after the merger is as- 
sumed to be small (Equation 140 1) . so that the change in spin 
after the merger due to accretion is negligible. 

In the left-hand panel, the mass of the more massive 
black hole is allowed any value in the range 1-4 x larger 
than the smaller black hole with a uniform distribution. 
The distribution therefore represents the final product of 



isotropic single-event major mergers, the distribution is cen- 
tred around (a) =0.58 and has a full width zero intensity of 
about 0.3. The right-hand panel represents equal-mass merg- 
ers only (1:1 mass ratio), it is centred around {a)=0.69 and 
is very narrow, with a full width zero intensity of about 0.05. 
The narrow distribution is partly due to the narrow initial 
spins, in the range 0-0.1. A broader distribution of initial 
spins would also lead to a broader distribution of final spins 
yet the mean would remain unchanged. 

We have carried out a second simple simulation, shown 
in Figure [131 Here, 10® black holes have each suffered a ran- 
dom number of mergers. All initial SMBH spins are drawn 
from a flat distribution between and 0.1, reflecting the 
low spin values expected from chaotic accretion. After the 
merger, they are assumed to accrete a small enough fraction 
of their own mass that the spin is not signiflcantly changed 
(the condition in Equation 1401 is not met). The exact num- 
ber of mergers each black hole will undergo is determined 
by drawing from a Poisson distribution with mean 0.7, re- 
flecting the mean numbe r of mergers the hos t galaxy has 
experienced since z ~ 1 l|Robaina et al.ll2010l ). The steUar 
masses of the bulges of the host galaxies are drawn from a 
Schechter function, and the masses of the merging black hole 
are assumed to be a flxed fraction of th is mass (1.4xlO~^ 
of the stellar mass, iHaring fc RbJ [20041). The final spin i s 
estimated using the formula given byl RezzoUa et al.[ l|2008l) . 
where we have randomised the orbital parameters (the orien- 
tations of each angular momentum component with respect 
to each other) as well as the initial black hole spins. 

In this simulation, 50% of the SMBHs do not undergo 
any mergers, and they retain the small initial spin (the fiat 
component around a=0-0.1). A significant fraction undergo 
1 or 2 mergers (35% and 12% respectively, peaking at 0.65 
and 0.8), and a small fraction undergo three or more mergers 
(3%, causing the high-spin tail). The distribution shows a 
sharp peak around a= 0.65, but the expectation spin is only 
(a) =0.35 (vertical dashed line in Figure [13[) . in very good 
agreement with our results from the best fitting set. This 
low mean value is due to the large fraction of black holes 
that have not undergone any mergers and have kept their 
original spins (0-0.1). The mean from this simple simulation 
is slightly lower than our estimate from fitting the radio 
LF of HEGs and LEGs, but the distribution is qualitatively 
similar. 

As discussed by IRobaina et al.l (|2010t ) , the mass limit 
of their survey requires both progenitor galaxies to have 
masses ^ 5 x 10^" M©. Hence major mergers involving one 
progenitor with a mass < 5 x 10^" M© are missed out, and 
their estimate is a lower limit on the mean number of ma- 
jor mergers experienced by galaxies with mass > 10" Mq. 
Indeed, studying the projected a xial distribution of galax- 
ies more massive than 10^^ Mm . [van der Wei et al.l (|2009l ) 
have argued that they must have all undergone a major 
merger. In addition, these mergers are expected to have 
been relatively gas-poor, otherwise the p rojected axial dis- 
tributions would be mo re similar to discs (|Naab et al.[[200^ : 
Ivan der Wei et a"i][2009l ). Hence the mean of 0.7 mergers per 
galaxy, on average, might be an underestimate. 

However, Figure [T3| does predict an approximately bi- 
modal distribution, qualitatively similar to what we infer at 
z=0 (see Section [5| and in particular Figure [6]), with rea- 
sonably similar expectation values. Our results suggest the 
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high-spin end to peak around a~0.80-0.99, rather than 0.7, 
but otherwise the distribution Figure[T3]is similar to the dis- 
tributions inferred for the LEGs. We do note that achieving 
the high spins we infer in our bimodal distributions, a~l, 
remains very difficult. 

The galaxy mergers driving the spin-up do not need to 
be completely 'dry' fe.g. Ivan Dokkurnll2005l ): provided that 
after the merger, the condition in Equation [40] is not met, 
the spin of the SMBHs will be predominantly determined 
by the coalescence of the two SMBHs as no significant spin- 
down will occur. Given a black hole mass of 10^ Mq needs 
to accrete a mass ~ 10* M© to spin down, and assuming a 
fraction 1.4x 10~'^ of th e galaxy's gas makes it to the black 
hole (e.g. lHaring fc Rixii2004 ). spinning down the most mas- 
sive black holes requires ~10^^ Mq of molecular gas. This is 
a very high gas mass, much higher than found in low-redshift 
elliptical galaxies (Lees et al. 1991; Young 2002). 

Another key issue is that our estimate of {a{z)) is 
dominated by black holes with masses ^10* M0, which 
reside in bulges with stellar masses ^W^^ M©, and these 
have generally elliptical morphologies. Thus, we have indi- 
rectly selected galaxies with elliptical morphologies, which 
is a proxy for selecting galaxies that have undergone at 
least one major merger. It is the fact that elliptical galax- 
ies have undergone major mergers that allows their central 
SMBH t o have high spin values. This point was noted by 
IScheueij (iJ19 92) and related ar g umen ts have also been pu t 
for ward bv IWilson fc Colbert! (jiggsl l. ISikora et all (|2007l ') 
and lVolonteri et al.l ^200^ . 

Figure [TJ] is only a simple simulation, and detailed 
studies of the cosmic evolution of SMBH spin s from merg 
ers can be found , for e x ample, in ,Vplpnteri et al.l (|200' 



iBerti fc Volonteril (|2008|'1. iLagos et al 



tions carried out by iFanidakis et al 



The Simula- 
incorporate a 

semi-analytic model of the accretion onto SMBHs in the 
centres of galaxies, and in this simulation mergers of galax- 
ies also cause mergers of the central SMBHs. The simula- 
tion tracks down the spin of the central SMBHs, and as- 
sumes chaotic accretion spins the black holes down. Fig- 
ure [14] shows the resulting spin history for black holes with 
m.^ 10* M0 (N. Fanidakis, 2010, priv. comm.). 

The behaviour seen in Figure [T^ is almost identical to 
our Figure [TT] This similarity extends also to the standard 
deviations, reflecting the broad distribution of SMBH spins 
at each redshift. The value of {a{z)) from the simulation is 
also in good quantitative agreement with our values at ZfZ,l- 

This line of argument suggests that their lower- 
luminosity AGN in spiral galaxies wi ll have SMBHs with 
low spin values (|King fc Pringlell2007l ). At first sight, this 
would appear in contradiction with the spin estimate for the 
nearby AGN MCG-6-30-1 a =0.98, from a broad iron line 
ijBrenneman fc Revnold3l2006l') . However, the h ost galaxy of 
this source is actually an SO (jFerruit et al.ll2000l ). so the high 
spin could have been acquired through merging. In addition, 
the caveats discussed in Section [7.21 should be kept in mind. 

Major mergers are also expected to happen at high red- 



" MCG-6-30-15 contains a SMBH with mass ~2xl0'3 Mq, ac- 
creting in the analogue to the high or very high state of GBHs. 
In fact, it is probably accreting close to the Eddingto n ratio, with 
A>0.8 jMcHardY et al.ll2005l : | Beckmann et al1l200^ . 
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Figure 14. Evolution of the mean spin as a function of red- 
shift for SMBHs with m, ^ 10* Mq, from the simulation of 
IFanidakis et al. The squares mark the mean value of 

the spin in each redshift bin, with the vertical bars representing 
the standard deviation, rather than error b ars. This incorp orates 
mergers of SMBH, using the prescription of lRezzoIIa et al. ( 2008]) 
as well as a prescription for 'chaotic accretion' llKinET ct al. 200^. 
The behaviour of the spin in the simulation is extremely similar 
to our inference of the evolution from observations (Figure [TTJ. 



shift {z ^ 1), but higher gas masses in the host galaxies 
imply higher accretion rates which will lead to a spin down 
(Equation 140 1 will be satisfied). Hence, while the mechanism 
to spin SMBHs up (major merging) is occuring both at high 
and low redshifts, the increase in spin at low redshifts can 
be attributed to a gradual disappearance of the 'braking' 
mechanism (chaotic accretion of a significant mass). 

If high-redshift quasars have a low mean spin due 
to chaotic accretion, it implies they have had time to 
spin down. The timescale for spin down is related to the 
timescale for alignment (~ lO'^ yr for a~l. Equation 15 of 
iLodato fc Pringlel [200^ ) since after a few alignments and 
anti-alignments the spin will have decreased to ~0.1. AGNs 
with a higher spin value are expected to have a longer 
alignment timescale (|Lodato fc Pringi3l2006l ). so jets from 
high-spin SMBHs will be more stable than those of SMBHs 
with low spin, implying powerful observed FRIIs c ould live 
for a few xlO^ yr (indeed the ag es estimated by Scheuei] 
1 19951 : IWillott et ai] Il999l : iBlundell fc Rawlinga ,2000l ) and 
will be rapidly spinning. Given that the most powerf ul radio- 
loud quasars typical jet advance speeds are~ 0.05c l|Scheueij 
1995), and assuming a maximum age of 5x lO'^ yr, the largest 
jets are expected to be ~8 Mpc at the most, in exc ellent 
agreement with observations (e.g. iBlundell et al.lll999l ). 

Our bimodal distribution of spins suggests that in the 
local Universe, the SMBHs with mass k, 10* Mq will be 
able to produce jets with only a modest rate of accretion 
(A^l). It is therefore no surprise that approximately 30% 
of local massive ellipticals have jets w ith moderate radio 
luminosities (e.g. lMauch fc Sadlerll2007h . 

Our estimate of {a{z)) is not in disagreement with the 
high values of a~0.9 inferred for the SMBHs of low-redshift 
elliptical galaxies (see N07). The authors inferred the spin 
from the correlation between Bondi accretion rate and jet 
power for nearby elliptical galaxies found by lAllen et al.l 
^20061 ), and using an ADAF model for the production of 
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jets around rotating SMBHs. The sources in the sample of 
lAllen et al.l (|2006) were selected to show jet-induced cavi- 
ties in the hot gas imaged using X-rays. Hence, the sample 
is biased towards elliptical galaxies with jets, which in our 
picture and that of N07 corresponds to high spin. Although 
our mean spin in the local Universe is typically lower than 
0.9, approximately half of the SMBHs accreting as ADAFs 
are predicted to have near-maximal spins (Section [4]), hence 
there is no conflict between our results and those of N07. 

Finally, we speculate how our inferred cosmic spin his- 
tory can be explained in the context of galaxy groups and 
clusters. We have argued that SMBH major mergers caused 
by galaxy major mergers will spin up the black holes. For 
major mergers to occur efficiently, there must be a high den- 
sity of massive galaxies, yet the translational velocities of 
the galaxies must not exceed the escape velocities of the 
galaxies themselves, which correspond to ~200-400 km s~^ 
for massive galaxies hostin g black holes of mass ^ 10* Mo 
l|Ferrarese fc Merrittllioool '). 

The ideal environment for mergers therefore corre- 
sponds to groups or to the outskirts of low-redshift clus- 
ters of galaxies. At z~ 2, the clusters of galaxies have not 
yet virialised fully, and the characteristic velocities of galax- 
ies in these high-redshift protoclusters will not exceed this 
velocity significantly. Hence, mergers are expected to be effi- 
cient, and the mean spin of SMBHs will be initially high due 
to mergers. However, at 2 large reservoirs of gas will also 
be present, allowing the SMBHs to accrete large masses. 

The most massive SMBHs will reside at the centres of 
these protoclusters: with black hole masses ~ 10^ M©, A > 
0.1 and initially o~ 0.6 due to mergers, these central SMBHs 
will produce powerful FRIT jets. The injection of energy from 
these FRII jets into the intracluster medium will prevent 
further cooling of the gas cutting down the supply of gas for 
future star-forma tion and SMBH growth in galaxies in the 
proto cluster (e.g. iBinnev fc Taboilll995l : iRawlings &: Jarvi3 
|2004| ). 

The SMBHs in the central dominant galaxies might 
well retain their high spin from the high-redshift epoch, 
z ~ 2, provided that the total mass accreted after their 
last major merger is small compared to the black hole mass, 
rn,~ 10^ Mq. Indeed the little evolution observed amongst 
the high mass end of the stellar mass function, with stellar 
mass >3 x 10^^ Mq, suggests little merging activity b elow 
2~1 for the most massive galaxies (|Banerii et al.|[2010f ). 

The bulk of the accretion is, however, in the form of 
radio-quiet quasars (m.~ 10* M©, A ^ 0.1 and a~0), whose 
spin at high redshift is low due to prolonged chaotic accre- 
tion. Presumably these objects have not yet undergone a 
major merger or have been able to accrete their own mass, 
Am,/m. ^ 1, to spin themselves down to a~0. The number 
of such objects is large compared to the number of ~ 10^ M© 
SMBHs in future central dominant galaxies, which brings 
the mean spin down. 

At 2: < 1 mergers of massive galaxies will still occur, 
particularly in environments with velocities ~300 km s~^, 
corresponding to group environments and the outskirts of 
clusters. These mergers will spin up the SMBHs. The gas 
in the host galaxies will have been stripped by ram pres- 
sure, while the intracluster gas is hot and will not resupply 
the galaxies. With little gas in the host galaxy, the braking 



mechanism, chaotic accretion of a significant mass, disap- 
pears. 

Future episodic accretion onto SMBHs with m,— 10* — 
lO" M0 will be modest and will not affect the spin. It will 
lead to intermitent FRI-ADAF activity, rather than to radio- 
quiet quasars. Occasionally a major accretion episode might 
occur, creating sources such as Cygnus A, but such episodes 
will be very rare. 
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